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ABSTRACT 


The Department of Defense has a continuing need for satellite communications to 
satisfy the demand for information exchange for strategic, operational, and tactical 
warfighters. There is currently a Deputy Undersecretary of Defense for Space (DUSD 
(Space)) transition planning effort to develop a satellite communications architecture for 
the 2007-2010 time frame. During this time all three current satellite systems; UFO, 
DSCS, and MILSTAR, are expected to degrade rapidly. As part of the DUSD - Space 
effort the U.S. Navy was tasked to form a Mobile Users Study to establish a framework 
for completing the detailed requirements and engineering work needed to develop the 
UHF/Mobile User transition plan. Then, as part of the Navy effort the Naval Postgraduate 
School's Astronautical Engineering class SE-61 under Professor Brij Agrawal's guidance 
designed a proposed medium Earth orbit communications satellite. This thesis is a design 
of the Attitude Dynamics and Control Subsystem for the subject medium Earth orbit 
MUS communications satellite. The thesis describes and explores the five major steps in 
designing an Attitude, Dynamics and Control Subsystem and focuses on key ADCS 


related areas that are peculiar to a MEO satellite as compared to a GEO satellite. 


CHOOL 
Of 


This Page Intentionally Left Blank 


Vi 


TABLE OF CONTENTS 


PP Tis) We) 0) MIO) rere eer scare ee etek oc. sce cscs csecesvscvcccecsscees 1 
MMII NAGS) EIN cs a erasers dcasssdcecenesscccccccccteescceese 1 

1. Mission Overview and Requirements...................0006. 1 

2. Management Responsibilities. .........0.......cccceeee eee eees + 

BPO ONCE PL Ol OMENAMONS: 6-2... ....c. ces cecececscecesesdeeusies 7 
PUMEMIACCCIAIE WCSCIIPUIOM a. +... ...00..0000ccsvieescodceccacsanes 9 


B. ATTITUDE DYNAMICS AND CONTROL 


SUBSYSTEM - CONCEPTS AND DESIGN PROCESS........ 16 
II. DEFINE CONTROL MODES AND SYSTEM 
He ony elem ine @NONC EVI INOS ot tiesec ste sie ees oc. sccscssveccesrenccseetaceresss 23 
free iN COON TROL MODES. ooo cc.00c..0ccseccssceteccesceens cess 23 
PPO IMIEMIISCT MOM VIOGC. 55.0 .2..5.c00t rere eee eee os ede e tet cease 25 
PP NCCUMSIUIOM WIOUE ee. isc ccdvseecscsiss.sccsdccetcesaserteees 23 
3. Sun Acquisition Mode / Sun Hold Mode.................... 23 
ame ATU Me NCCUUSILIOMIVIOUC.. 22.0.6. ..2..ccceccs cece tite sssss0ess 23 
Seer OU AlyON-OtAIOM WIOGC.. 252, -..0.5. 00 cceserdeecae++ sue 24 
GWM IG OC Mere ee he cok ve bovsu sls seuvddierss cesses 24 
Vs MOLONEY |) eS hy (076 | ae 24 
8. Orbit Maintenance / Orbit Transition Mode................ 24 
PSV eis Vie INO UTR EIMIIIN Docc. sce. cseceecccccsccssscseees 24 
Il. FUNCTIONAL DESCRIPTION AND TYPE OF SPACECRAFT 
CONTROL BY ATTITUDE CONTROL MODE......................065 oy 
A. ADCS FUNCTIONAL DESCRIPTION AND 
@ONMONE I OW ERWIE W o.oo ccecciceccccrccsccscersccesccves a 
NPM One TTT arenes sacs essed i dacscccescscees 29 
PRBS ACCOM MIC CIS: Regatta ele ancagcticec.. .sd0giidweeseee veces 2 


Vil 


3. Magnetic Torque Rodsay= Bs ee 29 
4, Star Sensor .::.......2:02e0. 015... 2. 1 eo 
5. Global Positioning System2--......5-- eee 30 
6. Inertial Reference Unit... .2... 3... 30 
7. Three Axis Mapnetemetcrgge,... ee 30 
8. Central Processor:......-2.. -34cucese ee 31 
9. Sun Sensor ..2.0..0.524.5. 9c oe eee ee 31 
10;.-Earth Sensors. 25.2 ...05. ... eee Bl 
B. TYPE OF SPACECRAFT CONTROL BY 
ATTITUDE. CONTROL MODE. 31 
lL. Orbit Insertion Mode........ 222-252 eee ee eee 31 
2. Acquisition Mode.. 2. 23a... -o-e eee ee 32 
3. Sun Acquisition Mode / Sun Hold Mode.................... 32 
4. Earth Acquisition Modes =...1-2---s neni ee 32 
5. Normal, On-Station Mode........7.22-.....245.5 eee 33 
6. Slew Mode. s:cccgyseseescs- vee eee 33 
7. Contingency / Safe Mode-..7.5eeeee eee 33 
8. Orbit Maintenance / Orbit Transition Mode................ 33 
C. THE MOTION OF THE SATELLITE AND THE 
SOLAR ARRAYS.........2.6 02222002 0c0c0ss cesses os ooo eee et ee 36 
IV. QUANTIFY DISTURBANCE ENVIRONMENT...............--220000: 43 
A. EXTERNAL DISTURBANCE TORQUEBG.................000000+ 43 
1. Solar Pressure Torques.)...jc94250 ee 44 
2. Gravity Gradient Torques.....29%:....... eee 45 
3. Aerodynamic Torques:.......4.5...0...-.. eee 46 
4. Disturbance Torques Due to the 
Geomagnetic Field:....2.05.25.4. ee 47 
B. MATLAB SIMULATION OF EXTERNAL 
DISTURBANCE TORQUES........:-2.0...2 eee eee 50 


Vill 


Pees Olatmencostite LOLGUES...¢............c.cccsceccecescseccceces 50 


Pmmeiavity Citadient POrgues....25............0c....scseceeecces 54 
Bee eromymanme i ressure UORQUES..........5.......0e0cceeceees 56 
4. Combined External Disturbance Torques................... 60 
C. INTERNAL DISTURBANCE TORQUEG..............000000000.. 63 


V. ATTITUDE DYNAMICS AND CONTROL SUBSYSTEM 
HARDWARE SELECTION AND DESCRIPTION............00.0000004. 65 
A. SIZING THE ADCS REACTION WHEELS. .................005. 65 
1. Torque and Angular Momentum Requirements 
eles oan? OUCeMIN Oa. .clees..5.. 24 Mee ss 65 


2. Torque and Angular Momentum Requirements 


Detar extemal MIStUrbANCeS,. 22h. iscsas geese. cese ee see 69 

B. SIZING THE MAGNETIC TORQUERSG..........0 00.0.0. 02 cece 69 
Geo MieheseNsOR SELECTION ISSUES.................0....04: 70 
Deis DWAR E ob Li Gane Nien... ....0RRaes.....ceccee ccs. ee 70 
lemspacecrait Control Processor..............0diis...essauudess 70 

2. Magnetic Acquisition and De-spin System................. 13 

3. Earth Sensor / Sun Sensor Assembly.....................06. 74 

MRCOG RVG Wines Mimneae ny aman foci eet hn, «'c 2 AEs vss: 74 

be Pheorno1 Operation. caer een 74 

4. The Reaction Wheel Assembly.......................00 cee 76 

A NPDICAHON’..... «..-. > ne NED, p NTI 5 76 

be Desion Description, yeameitt rs... 5.0otee eos. cs ascce 7] 

5. the satelite Inertial Reference System............:200200:- 79 

Glemthed © RIES C Mea tye ahs oiSate- cise y 4's «cs tanec ey Ceanaawe no 

b. Miniature Inertial Measurement Unit............... 79 

pM MMe Stal RACKED. a2. on, o).22 252-00 s.Csedsa cet eide on 83 
teGlobal Positioning Systemes ee 87 
ROU) SCNSOU MMMM ss iiissecce 2Ueu acess cease ooo OC 89 


b. Theory of Operation... 499303 89 
VI. YAW CONTROLLING ree 93 
A. YAW DISTURBANCES 23 repese rcs 0 sere eee ee 93 
B. SUN AND NADIR POINTING GEOMETRY.................... 93 
C. SUN AND NADIR YAW POINTING STRATEGY............. 94 
D. SUN AND NADIR ARRAY CONTROL STRATEGY.......... 96 
E. NORMAL, ON-STATION CONTROL FUNCTIONAL 
BLOCK. DIAGRAM eecce ees ee ee ee e 97 
F, ATTITUDE DETERMINATION SOFTWARE 
FUNCTIONAL OVERVIEW 2. ee 98 
G. SUN ANGLE ESTIMATION (B AND 8)...................00 0000 100 
H. YAW DATA PROCESSING Wty ee 101 


VU. SIMULATION OF THE ATTITUDE DYNAMICS AND 
CONTROL SUBSYS17 EM 222. ... . - Scenarios. .0 cena 
A. PROBLEM STATEMEN Dre recs rete ete err 
B.. THREE-AXIS STABILIZATION Ses. eee ee 
1. Small Angle Approximation.............. 2... ccccec eee eeee ens 
2. Large Angle Requirements 2227-29... .... eee ee 
C. THREE-AXIS REACTION WHEEL SYSTEM................... 
1 small Angle Approximation:25...5.2. eee... 
2. Yaw steening RequirementsS...7--2-..- eee 
D. ADCS CONTROL AND PERFORMANCE ANALYSIS - 
YAW STEERING. 2. ee. 


VIL, CONCLUSION vioicccc ssc See iets eee 
A. SUMMARY OP DESIGN PROCESS teas eee eee 
B. SUMMARY OF SPACECRAFT ADCS DESIGN................ 


forte NESTE Np OCOD YIM. 0 y5ce ccs essa cece ce secccccu ces sccccccscesscseeess 143 


EME UDP eee DO NIK pose se eisiglee sae ceccscccvescccessssccccsersececaacs 145 
APPENDIX C - MATLAB CODE FOR YAW STEERING.................. 149 
APPENDIX D - MATLAB CODE TO COMPUTE DCM .................... 155 
APPENDIX E - SIMULATION OF THE TOTAL EXTERNAL 
POMS ANG Sree. cc coe s cede cee scctscunddebesecccscans 157 
APPENDIX F - SMALL ANGLE SIMULATION ..........00.............000 169 
APPENDIX G - LARGE ANGLE SIMULATION.................. ccc cece ens Ns 
Pree) bette FIN INCE Sees esc cccs coe ccvecdicesecsescccccceedeabidvecersecere 181 
wee eI IIS UL TOIN LISI... eee eiitaee see cece eset lees cee ceee 185 


XI 


This Page Intentionally Left Blank 


Xil 


I. INTRODUCTION 


A. BACKGROUND 

The MEO UHF Satellite Constellation for Robust Assured Telecommunications 
(MUSCRAT) communications system design described in this chapter was created as a 
project for AA4871, Spacecraft Design II at the Naval Postgraduate School. The project 
provided experience in the design of a complex satellite at both the system and subsystem 
levels. Each student in the project team was assigned responsibility for a major subsystem 
or design support function. Design and analysis of each subsystem was performed over a 
period of eleven weeks using facilities available at the Naval Postgraduate School. 
Additional support was obtained from government and commercial sources as 
documented in this report. [Ref. 40] 

The MUSCRAT project was sponsored by the Naval Space Command. The goal 
of this project was to design a spacecraft bus and payload to meet the essential 
requirements for a Medium Earth Orbit (MEO) satellite communications system which 
will replace the aging UHF Follow-on (UFO) Constellation. Supporting trade studies, 
analyses, and options for increasing spacecraft capability are presented along with 


associated impacts on spacecraft mass, power, cost and schedule. [Ref. 40] 


1. Mission Overview and Requirements 


Advances in information technology place us at the point of a revolution in military 
operations — one in which the ability to collect, process, and disseminate information 
causes new thinking in doctrine, strategy, tactics, and procedures. Information and the 
dominant battlespace knowledge that it imparts allows a smaller, mobile joint force 
structure to be increasingly lethal and survivable in the battlespace. Success in modern 
warfare depends heavily upon achieving information dominance and that dominance 
depends upon the connectivity provided by communications systems. [Ref. 40] 

There is a compelling need for space-based communications assets to help meet 


the growing information needs of the advanced warfighting and supporting systems the 


Department of Defense (DoD) ts investing in for today and the future. Although land, sea, 
and air-based communications systems are vital to the warfighter, they cannot meet all of 
the warfighters’ information needs. Space based communication assets, due to their instant 
accessibility, survivability, coverage, flexibility, and ability to support mobile subscribers, 
are in some cases the only practical means to support mobile warfighter platforms, such 
as ships, submarines, aircraft, and land combat vehicles. As a result, the demands on 
space-based communications are growing. Recent DoD studies and industry assessments 
project at least a three to five-fold growth in warfighter information requirements that 
must be satisfied by SATCOM systems in the next decade. [Ref. 40] 

Under the direction of the Deputy Under Secretary of Defense for Space 
(DUSD(Space)), a Mobile Users Study (MUS) was commissioned on 19 February 1997 
to define the requirements, architecture, and cost for a replacement to UFO. The Naval 
Space Command was tasked as the lead for the Navy-directed MUS Requirements 
Working Integrated Product Team (WIPT). This product team, consisting of 52 members 
from all military staffs, the Jomt Chiefs of Staff (JCS-J65), the National Reconnaissance 
Office (NRO), and the Defense Information Systems Agency (DISA), reviewed over 100 
general requirements for low data rate netted communications and narrowed them into a 


basic set of 8 high level needs: 


a. Assured Access 


Assured access to SATCOM services is the most fundamental SATCOM 
need of the warfighter. The SATCOM study process confirmed that deployed and mobile 
warfighters are largely dependent on SATCOM to satisfy their most critical beyond-line- 
of-sight information transfer needs. Unified Combatant Commanders have repeatedly 
requested the ability to access SATCOM on demand and to control the resources 
apportioned to them by the Joint Chiefs of Staff (JCS). A warfighter’s access to 
SATCOM support must therefore, be available on-demand when and where needed for 
the duration of the mission. [Ref. 40] 


b. Netted Communications 


A communications topology should be provided that allows multiple users 
to transmit (one at a time) and receive (all users less the one transmitting) with others on a 
common data stream. Netted communications on demand exists when flexible and 
reconfigurable nets enable the net control to bring on new members to the net without 


difficulty. [Ref 40] 


c. Communications on the Move 

This capability shall provide the ability of the warfighter to move, and talk, 
at the same time. This capability enables the user to have voice and data communications 
in any natural environment to include a double canopy / jungle, in rain at rain rate mode 


“H”, and in a sea environment. [Ref. 40] 


d. Jointly Interoperable Communications 
Current and emerging warfighting doctrine require joint operations to 
include inter-service compatibility and operating with allies. Joint service and allied 


operations require jointly interoperable communications. [Ref. 40] 


e. Worldwide Coverage 


Worldwide coverage is defined as24 hours a day communications service 
from 65 degrees N to 65 degrees S_ latitude without gaps in geographical coverage. [Ref. 
40] 


f. Point to Point Communications 

The system shall provide a communications topology that allows a single 
terminal to connect to another single terminal. Examples are communications link between 
a commander and a single subordinate element, a control platform or a remote device, or 


even a normal telephone connection. [Ref. 40] 


g. Broadcast Communications 

The system shall provide a communications topology that allows a single 
user to transmit while multiple users receive. A “one to many” broadcast capability is 
particularly valuable for the transmission of intelligence products or any information that 


needs to be distributed from a single source to multiple users. [Ref. 40] 


h. Polar Coverage 


Polar coverage is defined as 24 hours a day communications service for the 
area of the Earth’s surface that is above 65° N and below 65° S latitude. Communications 
service coverage primarily directed at the northern polar regions is important for reasons 


of national security. [Ref. 40] 


2. Management Responsibilities 

Management functions and responsibilities for MUSCRAT will mirror those 
assigned to current Military Satellite Communications (MILSATCOM) systems as set 
forth in CJCS MOP-37. Figure 1-1 illustrates the organizational structure for 
MUSCRAT. Table 1-1 outlines the responsibilities of each of the organizations involved. 
[Ref. 40] 
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Figure 1-1: Organization Structure for MUSCRAT 






Joint Staff 1. Approve requirements as part of the Integrated Communications Data 
Base (ICDB) process 
2. Apportion satellite resources to CINCs 
3. Adjudicate resource apportionment conflicts not resolved at lower levels 
4. Approve the Operational Requirements Document (ORD) 


Defense Information 1. Administer the ICDB 
Systems Agency (DISA) 


Unified Commanders-in- Review, validate, and forward requirements to the Joint Staff 
Chief (CINCs) Collect, validate, and prioritize satellite access requests from their 
components 

Manage satellite resources apportioned to them by the Joint Staff 
Provide management oversight of communications payload control at 
the direction of the Joint Staff and in coordination with the System 
Manager 
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Commander-in-Chief, 
United States Space 
Command 
(USCINCSPACE) 


Chief of Naval Operations | 1. Carry out the Department of the Navy’s responsibilities as the System 
(CNO N63) Manager for MUSCRAT 
Commander, Naval Space | 1. COMNAVSPACECOM is MUSCRAT’s system operational manager 


Command 2. Develop and coordinate procedures for payload control 

(COMNAVSPACECOM) | 3. Manage communications payload configuration to support CINC 
requirements 

4. Review requirements and perform analyses to determine feasibility of 
satisfaction via the MUSCRAT system 

5. Recommend apportionments to the Joint Staff 

6. Provide system use and status information to the Joint Staff, DISA, and 

USCINCSPACE 

Coordinate MUSCRAT payload health and welfare with the AFSCN 

ground stations and perform payload control functions 

2. Perform anomaly detection and resolution in conjunction with the 
AFSCN ground stations 

3. Actas controlling authority for the transmission security (TRANSEC) 

crypto key system, and direct uploading of the new key segment 

variables to the onboard TRANSEC device as required 

Coordinate with NAVSOC concerning satellite health/welfare and 

anomaly detection/resolution 
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Air Force Satellite l. 
Control Network 
(AFSCN) 
Naval Computer and 
Telecommunications Area 
Master Stations 
(NCTAMS) and Joint 
Fleet Telecommunication 
Operations Centers 
(JFTOCs) 


Users 


























1. Coordinate with NAVSOC concerning satellite payload control 
2. Provide connectivity to the Defense Switching Network (DSN) and 
Public Switching Telephone Network (PSTN) 






















1. Obtain approval of their requirements through the ICDB process 
2. Adhere to the network/terminal configurations and precedences assigned 
by the CINC communications manager 

Ensure that payload resources are used as assigned 


Table 1-1: Responsibilities of Organizations Involved with MUSCRAT 
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3. Concept of Operations 

MUSCRAT will be launched out of Cape Canaveral, Florida onboard either a 
Delta III or one of the Atlas family of launch vehicles (IIAS, ITAR, ITARS). Operational 
checkout is expected to last for approximately one month, after which time, the satellite 
will be placed in a non-operational mode until the entire constellation has been placed in 


orbit. The MUSCRAT system will then be activated. 





Figure 1-2: MUSCRAT Concept of Operations 


Figure 1-2 illustrates MUSCRAT’s concept of operations. The system uses 
current Air Force Satellite Control Network (AFSCN) stations for controlling satellite bus 
operations, as well as for monitoring standard satellite telemetry parameters. Control of 
payload operations will be accomplished through the current Naval Computer and 
Telecommunications Area Master Stations (NCTAMS) = and_ Joint’ Fleet 
Telecommunications Operations Centers (JFTOCs) via an S/EHF link. The MUSCRAT 


design allows for connectivity of payload control commands when a satellite 1s not within 
line of sight of a NCTAMS/JFTOC via EHF satellite crosslinks (60/61 GHz). 
MUSCRAT provides direct UHF connectivity between users (uplink: 290- 

320 MHz, downlink: 243-270 MHz). If one of the users are non-mobile and is using the 
Defense Switching Network (DSN) or the Public Switching Telephone Network (PSTN), 
gateways located at the NCTAMS/JFTOCs will be used. This greatly increases 
MUSCRAT?’s capacity and flexibility. The gateways connect to the satellite via an S/EHF 
link (downlink: 20 GHz, uplink: 44 GHz). The EHF satellite crosslinks will also be used 
for payload communications connectivity. Examples of these different modes are shown 


in the figure. 


4. Spacecraft Description 

MUSCRAT will provide an assured, netted, global (including the North and South 
Poles) UHF communications link between both mobile and fixed-site users. The 
MUSCRAT constellation consists of ten operational satellites divided into two orbital 
planes separated by 180° in longitude of ascending node. Their circular orbits have an 
altitude of 11,000 km and an inclination of 48 degrees. MUSCRAT can be launched on 
Atlas ITAR, Atlas I[AS, and Delta III launch vehicles. The constellation guarantees global 
coverage with at least one satellite within view of a user with a ten degree elevation angle. 
[Ref. 40] 

Figures 1-3 through 1-5 depict MUSCRAT’s deployed and stowed configurations. 
MUSCRAT is basically rectangular in shape, with a bus height of six meters and a base 
measuring two meters square. The panels for the main “box” are constructed of aluminum 
honeycomb and are supported with aluminum shear webs. MUSCRAT also uses an 


aluminum central cylinder for load support. [Ref. 40] 






Figure 1-3: MUSCRAT Deployed Configuration 
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Figure 1-4: MUSCRAT Deployed Configuration 
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Figure 1-5: MUSCRAT Stowed Configuration Inside Delta Fairing 





The bus has a modular design (Figure 1-6). The upper module contains the UHF 
components. The middle module houses the EHF downlink components, solar arrays, 
attitude control equipment, and control electronics. The lower module includes the EHF 
satellite crosslink components and batteries. The aluminum central cylinder runs the 
length of the bus, and encloses one pressurant and two propellant tanks (Figure 1-7). 


(Ref. 40] 


UPPER MODULE: 
UHF Payload Components 
Heat Pipes 


MIDDLE MODULE: 
Central Computer 
EHF Gateway Components 
ADCS Components 
EPS Solar Arrays and 
Control Electronics 


LOWER MODULE: 
EHF Satellite Crosslink Components 
EPS Battenes 





Figure 1-6: MUSCRAT's Modular Bus Design 


The most prominent feature of MUSCRAT is its 8.5 meter diameter reflector. It 
deploys on a graphite truss structure 5.3 meters above the bus. The ground plane, with a 
diameter of four meters, contains 19 crossed dipole antennas. Twenty-six high heat 
dissipating solid state power amplifiers (SSPA’s) are mounted on the north and south 
panels of the spacecraft. Seven of these are for redundancy. The remaining UHF 


components, except the data processing units, are located throughout the upper module. 











Figure 1-7: MUSCRAT's Central Cylinder 


MUSCRAT uses an EHF link to connect to gateway ground stations. These 
stations are used to link communications to the Defense Switching Network (DSN) and 
the Public Switching Telephone Network (PSTN). Two EHF antennas (plus one for 
redundancy) accomplish this connection. MUSCRAT also uses EHF satellite crosslinks to 
ensure timely netted communications. Four EHF 0.45 meter dish antennas are located on 
the lower corners of the bus. 

Control of spacecraft operations, with the exception of attitude determmation and 
control, is accomplished using a central data processor. It contains eight Motorola Power 
PC 603e microprocessors mounted on modified Radstone PCP2 conduction cooled 


integrated circuit cards. 





MUSCRAT is three-axis stabilized and yaw steered. Pointing requirements of 0.1° 
are met using three reaction wheels, with a fourth reaction wheel for redundancy. 
Positional and three axis attitude data are acquired using GPS and star trackers. The 
reaction wheels are desaturated using magnetic torque rods and 0.5 Ib thrusters as primary 
and secondary systems respectively. The attitude determination and control subsystem 
(ADCS) is controlled by its own redundant processor. Other supporting ADCS 
components include a magnetometer, two inertial reference units (IRU’s), sun trackers 
located on the solar arrays, and two earth/sun sensor assemblies used to initially orient the 
spacecraft once it gets to orbit. 

The telemetry, tracking, and control (TT&C) subsystem provides for both 
autonomous operations and direct ground control of MUSCRAT. The satellite has two 
patch TT&C omni-directional antennas for assured telemetry download and command 
reception. Onboard TT&C operations will be controlled by redundant spacecraft control 
electronics. TT&C ground functions will take place at Air Force Satellite Control 
Network (AFSCN) Space-Ground Link System (SGLS) stations. 

The electric power subsystem (EPS) is a dual bus, with a fully regulated bus 
voltage of 51 volts. MUSCRAT requires 8083 W of power in non-eclipse periods. Two 
GaAs solar arrays with a combined area of 30 m’ provide this power, and it is stored in 
four NiH, batteries, each having 27 cells. Shunt regulators are located on the solar array 
yokes. 

The propulsion subsystem is comprised of two monopropellant (hydrazine) tanks 
and a nitrogen pressurant tank housed inside MUSCRAT’s central cylinder. Propellant 
will be used for satellite re-phasing, orbital insertion corrections, as a backup for reaction 
wheel desaturation, and for satellite deorbit at the end of mission life. Twelve 0.5 lb 
thrusters are used for reaction wheel desaturation, and two 5.0 lb thrusters are used for 
gross orbit maneuvers. All thruster valves have dual series seat valves for redundancy 

MUSCRAT’s thermal subsystem provides the required temperature control 
throughout the satellite’s lifetime. The primarily passive system design uses heat pipes in 
the north and south panels to cool concentrated heat sources (i.e. SSPA’s) and minimize 


panel temperature gradients. These panels are radiators, and are covered with OSR. 


Multi-layer insulation (MLI) blankets are also used to maintain proper temperatures. 
Heaters are used to maintain the batteries and hydrazine system within required 
temperature limits. 

Table 1-2 is a summary of MUSCRAT’s mass budget. The total launch mass of 
2380.5 kg is broken down by subsystem, and includes a 15% mass margin, propellant, and 


the spacecraft adapter. 


110.1 


15% Margin 
Spacecraft Adapter 


Table 1-2: MUSCRAT Mass Summary 





Table 1-3 summarizes MUSCRAT’s power budget during eclipse and non-eclipse 
periods. It includes a 5% growth margin in sun, a 10% growth margin in eclipse, and an 


additional 10% solar array design margin. 
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* Not included in total power requirements due to extremely low average 
power requirements 


Table 1-3: MUSCRAT's Power Budget 


Table 1-4 lists MUSCRAT’s propellant budget. The total propellant mass of 
246.1 kg includes a 10% margin. 









Orbit Injection Correction 12.69 
Orbit Rephasing 72.54 


Mass (kg) 
Reaction Wheel Desaturation 17.28 






118.98 
10% Margin 24.61 
Total Propellant Mass 246.1 


Table 1-4: MUSCRAT's Propellant Budget 


B. ATTITUDE DYNAMICS AND CONTROL SUBSYSTEM - 
CONCEPTS AND DESIGN PROCESS 


The attitude determination and control subsystem (ADCS) stabilizes the vehicle 
and orients it in desired directions during the mission despite the external disturbance 
torques acting on it. This requires that the vehicle determine its attitude, using sensors, 
and control it, using actuators. The ADCS is often tightly coupled to other subsystems on 
board, especially the propulsion and navigation functions. [Ref. 6: p. 340] 

This paper begins by discussing several useful concepts and definitions, including 
mass properties, disturbance torques, angular momentum, and reference vectors. The 
mass properties of a spacecraft are key in determining the size of control and disturbance 
torques. And, a knowledge of the location of the center of mass or gravity (cg) as well as 
the elements of the inertia matrix is a requirement. As well as the moments and products 
of inertia about chosen reference axes. The direction of the principal axes, those axes for 
which the inertia matrix is diagonal and the products of inertia are zero, are also of 
interest. Finally, we need to know how these properties change with time, as fuel or other 
consumables are used, or as appendages are moved or deployed. [Ref. 6: p. 340] 

A body in space is subject to small but persistent disturbance torques (Le., 
10° *Nm) from a variety of sources. These torques are categorized as cyclic, varying in a 
sinusoidal manner during an orbit, or secular, accumulating with time, and not averaging 
out over an orbit. These torques would quickly reorient the vehicle unless resisted in 
some way. An ADCS system resists these torques either passively, by exploiting inherent 
inertia or magnetic properties to make the "disturbances" stabilizing and their effects 
tolerable, or actively, by sensing the resulting motion and applying corrective torques. 
[Ref. 6: p. 340] 

Angular momentum plays an important role in space, where torques typically are 
small and spacecraft are unconstrained. For a body initially at rest, an external torque will 
cause the body to angularly accelerate proportionally to the torque - resulting in an 
increasing angular velocity. Conversely, if the body is initially spinning about an axis 


perpendicular to the applied torque, then the body spin axis will precess, moving with a 


constant angular velocity proportional to the torque. Thus, spinning bodies act like 
gyroscopes, inherently resisting disturbance torques in two axes by responding with 
constant, rather than increasing, angular velocity. This property of spinning bodies, called 
gyroscopic stiffness, can be used to reduce the effect of small, cyclic disturbance torques. 
This is true whether the entire body spins or just a portion of it, such as a momentum 
wheel or spinning rotor. [Ref. 6: p. 340] 

Conservation of vehicle angular momentum requires that external torques change 
the system angular momentum. Thus, external disturbances must be corrected by external 
control torques (i.e., thrusters or magnetic torquers) or the resulting momentum buildup 
must be stored internally (i.e., by reaction wheels) without reorienting the vehicle beyond 
its allowable limits. The momentum build up due to secular disturbances ultimately must 
be reduced by applying compensating external control torques. [Ref. 6: p. 340] 

Often, in addition to rejecting disturbances, the ADCS must reorient the vehicle (in 
slew maneuvers) to re-point the payload, solar arrays, or antennas. These periodic re- 
poimting requirements may drive the design to larger actuators than would be required for 
disturbance rejection alone. [Ref. 6: p. 341] 

To orient the vehicle correctly, external references must be used to determine the 
vehicle's absolute attitude. These references include the Sun, the Earth's IR horizon, the 
local magnetic field direction, and the stars. In addition, inertial sensors (gyroscopes) also 
can be carried to provide a short-term attitude reference between external updates. 
External references (1.e., Sun angles) are usually measured as body-centered angular 
distances to a vector. Such measurements provide only two of the three independent 
parameters needed to specify the orientation of the spacecraft. This results in the need for 
multiple sensor types on board most spacecraft. [Ref. 6: p. 341] 

The Attitude Dynamics and Control Subsystem design process is summarized very 
nicely with Table 1-5. Much of the input information for this table comes from other 
subsystems and the overall concept of operations for the satellite. This can be seen in 


Figure 1-8. [Ref. 6: p. 344] 
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Table 1-5: Attitude Dynamics Control Sub-System Design Process 
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Figure 1-8: The Impact of Mission Requirements and Other Subsystems on the ADCS 


Subsystem. Direction of arrows show requirements flow from one subsystem to another. 


The first step in designing the Attitude Dynamics and Control Subsystem is to 
define the control modes. Each control mode is related to the different phases of the 
spacecraft life. The ADCS for this satellite will have eight control modes which are 
described in chapter 2. The sensor and actuator requirements for each control mode may 
vary drastically. As a result, the designer of the ADCS must ensure the requirements of 
each control mode are satisfied while at the same time producing as simple and reliable a 
subsystem as possible. 

Next, a thorough understanding of the system level requirements is required. 
These requirements usually originate from the other spacecraft subsystems and flow down 
to the ADCS designer. For example, the payload designer in this design specified the 
UHF antenna must maintain a pointing accuracy of 0.1°. Then, the systems engineer and 
structures designer determined the UHF antenna would be fixed to the spacecraft. This 
can then be translated into a spacecraft attitude control accuracy of 0.1° along the roll and 
pitch axes. 

Then, based on mission requirements and inputs from other subsystems, a decision 
has to be made on the type of spacecraft control to be used. For example, a spin stabilized 
spacecraft may be most suitable for some missions while a three-axis control momentum 
wheel system may be best for others. When the type of spacecraft control is determined 
an understanding of the type of sensors and actuators should be developed. This is not 
actually selecting the components, but instead performing trade studies for determining 
which types of sensor and actuator are best suited for the spacecraft mission. An example 
might be determining whether to use a star sensor or earth sensor. Chapter three 
discusses these issues and presents a functional description of the attitude dynamics and 
control subsystem. 

The next step in the design process is to quantify the disturbance environment. 
There are generally four external disturbances to consider: Gravity Gradient, Solar 
Pressure, Aerodynamic Pressure and Geomagnetic field disturbances. These disturbances 
should be estimated so that appropriate actuators and sensors can be selected to 
counteract their effects. Chapter four presents a thorough explanation on quantifying the 


external disturbance torques. 
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With information from the analysis discussed above a selection of ADCS 
components was performed. The components of the ADCS include sensors, actuators, 
and control processors. Chapter 5 presents this selection process and provides a detailed 
description of each component selected. 

The final portion of the design process is the definition of determination and 
control algorithms. This also includes simulation and modeling of the ADCS. This begins 
in chapter six which focuses exclusively on yaw controlling and control algorithms. 
Chapter 7 then focuses on simulating the motion of the spacecraft and presents a thorough 
explanation of the differences between a non yaw Steering satellite and a yaw steering one. 


The thesis concludes with several spacecraft attitude control simulations. 


21 





This Page Intentionally Left Blank 


Ze 


fl. DEFINE CONTROL MODES AND SYSTEM LEVEL 
REQUIREMENTS 


Attitude control modes are related to the different phases of the spacecraft life. 
Basically, the ADCS must be designed to meet the requirements of each mode. These 
ADCS requirements are tightly coupled to other spacecraft subsystems and the overall 
mission needs. Since the requirements change drastically with each mission phase, the 


ADCS designer may be challenged to built one system that meets all requirements. 


A. DEFINE CONTROL MODES 

1. Orbit Insertion Mode 

The orbit insertion is the period during and after boost while the spacecraft is 
brought to final orbit. Options in this control mode include no spacecraft control, simple 
spin stabilization of solid rocket motor, and full spacecraft control using liquid propulsion 


system. 


2. Acquisition Mode 


This is normally where the initial determination of attitude and stabilization of the 
vehicle occurs. But, attitude determination and stabilization may also need to occur due 


to recovery from power upsets and / or emergencies. 


3. Sun Acquisition Mode / Sun Hold Mode 


After attitude determination and stabilization is complete Sun acquisition takes 


place in order to ensure an adequate power supply. 


4. Earth Acquisition Mode 
A critical aspect of this mode is acquiring a nadir pointing payload orientation 
while simultaneously keeping the solar arrays tracking the sun. Information from the sun 


sensors and additional navigational aids allow the spacecraft to re-orient such that 


ZS 


simultaneous Nadir payload pointing and solar array sun tracking is attained. Much more 


will be discussed on this topic later. 


5. Normal, On-Station Mode 


The Normal, On-Station mode is used for the vast majority of the mission. As a 
result, the requirements for this mode will drive the system design. And, the majority of 


discussion will focus on this mode of operation. 


6. Slew Mode 


This mode occurs when a reorientation of the vehicle is required (i.e. the 


spacecraft executes a slew maneuver). 


7. Contingency / Safe Mode 
This mode of operation is used in emergencies for when the regular mode fails or 
is disabled. Also, this mode may be selected to meet unexpected requirements like power 


or thermal constraints; at the expense of normal operations. 


8. Orbit Maintenance / Orbit Transition Mode 


This mode occurs when the spacecraft needs to be repositioned in the orbit. This 
either could be due to an accumulation of error over several years; or due to a satellite 


failure which requires selected satellites to be repositioned in the orbit. 


B. SYSTEM LEVEL REQUIREMENTS 


Following are the system level requirements for the Attitude Dynamics and 
Control Subsystem of the MUS communications satellite. 

1. The ADCS must be fully operational without the aid of ground base control. 

2. The operational lifetime of the ADCS must meet or exceed ten years. 

3. The ADCS must maintain a UHF antenna beam pointing accuracy of less than 


or equal to 0.1°. 
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4. The ADCS must have the capability to provide continuous yaw steering and 
solar array tracking such that the solar array pointing accuracy is better than or equal to 
5: . 

5. The spacecraft must be fully operational during eclipse. 

6. The spacecraft must be stabilized and under control during all phases of the 
mission. 

7. The ADCS system should have system redundancy such that there are no single 
point failures. 

8. The ADCS should have a ground override capability. 

9. The spacecraft must have a cross-link antenna pointing accuracy of 0.07° which 
requires a spacecraft attitude knowledge requirement of 0.02° in roll, pitch, and yaw. 

10. The Attitude Dynamics and Control Subsystem must have a positional 


knowledge requirement of 350 meters. 


Zs 
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Hil. FUNCTIONAL DESCRIPTION AND TYPE OF SPACECRAFT 
CONTROL BY ATTITUDE CONTROL MODE 


A. ADCS FUNCTIONAL DESCRIPTION AND COMPONENT 
OVERVIEW 


A Zero Momentum three axis system was selected which can provide any 
combination of slew maneuvers in the pitch, roll and yaw directions. The Attitude 
Dynamics and Control Subsystem will maintain Nadir pointing accuracies to within 0.1°. 
At the same time it will continuously rotate in the yaw direction while the solar arrays are 
simultaneously rotating in order to keep a perpendicular orientation towards the sun. This 
will greatly decrease the complexity of the solar array control mechanisms. In a zero- 
momentum system, reaction wheels respond to disturbances on the vehicle. For example, 
a vehicle-pointing error creates a signal which speeds up the wheel, initially at zero. This 
torque corrects the vehicle and leaves the wheel spinning at low speed, until another 
pointing error speeds the wheel further or slows it down again. If the disturbance is cyclic 
during each orbit, the wheel may not approach saturation speed for several orbits. Secular 
disturbances, however, cause the wheel to drift toward saturation. We then must apply an 
external torque with a thruster, a magnetic torquer or both, to force the wheel speed back 
to zero. This process is called desaturation, momentum unloading, or momentum 
dumping and can be done automatically or by command from the ground. [Ref. 6: p. 348] 

The Zero Momentum Three Axis Control system chosen was one of three options 
seriously considered. The other two were Dual Spin Stabilization and a Zero Momentum 
Control Momentum Gyro. The Dual Spin Stabilization was dismissed because the antenna 
configuration did not fit well on the despun platform and the power requirements were 
greater than a Dual Spin Stabilization could provide. The Control Moment Gyro (CMG) 
would work for this system, but it was dismissed because it is heavier and would not be 
conducive to yaw steering. Yaw steering enables a decrease in mass and complexity of the 
solar array drive assembly. In addition the Zero Momentum Three Axis Stabilization 


system would be lighter. Also, the Control Moment Gyro is best suited for high torque 
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demands and fast slewing requirements. Neither of which we have. It should be pointed 


out however, that the three axis stabilization has more moving parts and could therefore 


be considered more complex. So, a legitimate argument could be made for either system. 


Figure 3-1 is a general block diagram of the Zero Momentum Attitude Dynamics 


Control Subsystem. [Ref. 6: p. 348, Ref. 21: p. 46] 
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Figure 3-1: Block Diagram of Attitude Determination and Control Subsystem 
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1. TT&C 


Telemetry, Tracking and Control will take place at the seven Air Force Satellite 
Control Network ground stations with one of the stations exercising overall control. 
Signals will be sent to the spacecraft as required where they will be processed in the 
Spacecraft's main computer. Those commands relating to Attitude Dynamics and Control 
will be sent to the ADCS central processor and then executed. In addition, data will be 
sent from the ADCS system to the ground stations in order to monitor the health, welfare 
and position of both the satellite and the ADCS. The TT&C subsystem will consist of 
dedicated antennas, a receiver, a transmitter, and a command and data handling system. 


[Ref. 15: p.13] 


2. Reaction Wheels 


The Zero Momentum reaction wheel system will counteract disturbance torque 
and provide the momentum to perform slew maneuvers as required. Reaction wheels are 
characterized by being light weight and capable of producing a continuous range of 
torque. This will work well for our satellite sce we will continuously need to track the 


sun in order to attain maximum capacity with our solar arrays. [Ref. 21: p. 47] 


3. Magnetic Torque Rods 

Magnetic Torque Rods will be used to efficiently dump angular momentum in the 
reaction wheels produced by secular disturbance torques. By providing electric current 
through the torque rod coil a magnetic dipole moment results. This moment causes a 
torque that tries to align the torque rod with the earth's magnetic field. When this torque 
is produced an equal and opposite torque will be produced by the Reaction Wheels; thus 


allowing a means to dump momentum. [Ref. 21: p.55] 
4. Star Sensor 


Based on the crosslink antenna pointing requirements of 0.02° some type of star 


sensing system will have to be used. In order to provide a continuous star sensing three 
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star sensors will be required. In the worst case scenario two of the three sensors may be 
obstructed by the sun and moon. If this occurs the Inertial Reference Unit will be able to 


provide the required data for the short period of time required. 


5. Global Positioning System 


A Global Positioning System (GPS) is required on the spacecraft in order to allow 
for an ability to accurately determine orbit position without the aid of a ground station. 


And, the GPS will be a back-up for attitude determination. [Ref. 21: p. 53] 


6. Inertial Reference Unit 


The Inertial Reference Unit will provide continuous angular velocity information. 
Then, with the knowledge of angular velocity, angular acceleration can be obtained using 
a pseudorate modulator. In the event of star tracker failure or temporary un-availability, 
the Inertial Reference Unit will also provide continuous attitude information for a limited 


period of time. [Ref. 21: p. 52] 


7. Three Axis Magnetometer 


Magnetometers are simple, reliable, lightweight sensors that measure both the 
direction and size of the Earth's magnetic field. When compared to the Earth's known 
field, their output helps establish the spacecraft's attitude. But their accuracy is not as 
good as that of star, sun or horizon references. The Earth's field can shift with tume and 
cannot be calculated precisely in the first place. To improve accuracy, data is often 
combined with the data from Sun or horizon sensors. When the vehicle is using the 
magnetic torquers and passes through magnetic field reversals during each orbit, the 
magnetometer will be used to control the polarity of the torquer output. And, the 
torquers will have to be turned off while the magnetometer is sampled to avoid corrupting 


the measurement. [Ref. 6: p. 361] 
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8. Central Processor 

The Attitude Dynamics and Control Subsystem Central Processor will be a 
commercial off the shelf product produced by Honeywell. The Central Processor is the 
brain of the ADCS. As can be seen from Figure 3-1, input data is fed into the central 
processor from the TT&C, Star Sensors, GPS, IRU, Three Axis Magnetometer and 
Reaction Wheels. Based on this input data the central processor will calculate the 
spacecraft's position, velocity, orientation, angular velocity and angular acceleration. 
Then, the central processor will calculate the required spacecraft orientation, angular 
velocity and angular acceleration profile to maintain Nadir earth pointing, an adequate sun 
pointing profile, momentum limits and orbit repositioning/transfers. Finally, appropriate 
commands will be sent to the reaction wheels, torque rods, thrusters and solar array drive 


mechanisms. 


9. Sun Sensor 
A sun sensor is a visible light detector which can measure one or two angles 
between their mounting base and incident sunlight. One sun sensor assembly will be 


mounted on or near each solar array assembly and will be used to track the suns rays. 


10. Earth Sensor 
An Earth Sensor is an infrared device that detects the contrast between the cold of 
deep space and the heat of the Earth's atmosphere. These sensors provide Earth relative 


information directly for earth pointing spacecraft, which simplifies on board processing. 


B. TYPE OF SPACECRAFT CONTROL BY ATTITUDE CONTROL 
MODE 
1. Orbit Insertion Mode 


The launch vehicles we have selected will allow for direct orbit insertion with no 


need for an apogee kick motor or a perigee kick motor. As a result, under normal 
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condition there will be no spacecraft control until after the spacecraft is in or very near its 


final orbit and has separated from the launch vehicle. 


2. Acquisition Mode 


When the spacecraft is separated from the launch vehicle the acquisition mode of 
operation will begin. At this time the spacecraft will normally have a spin about the z axis. 
Each of the two sun sensor assemblies have two perpendicular fan shaped fields of view, 


a, and a,. Each sensor generates an output pulse as the sun passes through the field of 


view Of the light sensitive element. Then, the time of arrival of the pulse is detected in the 
ADCS central processor. In addition, the earth sensor provides a pulse which is also 
detected by the central processor. The central processor is then able to take these 
measurements and compute the attitude of the spacecraft with respect to the sun and 


earth. 


3. Sun Acquisition Mode / Sun Hold Mode 

Using angular position, angular rate and data from the pitch sun sensor, the 
spacecraft is able to lock into and hold the sun in pitch. At this time the yaw and roll axes 
are under gyro reference control. Then, a slow yaw rotation occurs in order to acquire the 
sun with the yaw sun sensor. When both sensors have locked onto the sun, the sun hold 
mode automatically begins. Then, with the sun's position firmly established, the solar 
arrays are deployed and the remaining ADCS components are brought on line and 


thoroughly tested. 


4. Earth Acquisition Mode 


The earth acquisition mode occurs at either 0600 or 1800 local time. At this time 
the earth line and sun line relative to the spacecraft are perpendicular. With the yaw and 
pitch axis locked onto the sun a gentle roll maneuver takes place until the earth is acquired 
by the earth sensors. When an adequate Nadir pointing accuracy is attained by the earth 


sensors, both roll and pitch control transfer from the Gyros over to the earth sensors. 
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5. Normal, On-Station Mode 


The Normal/On-Station operating mode is by far the most significant of the 
Attitude Control Modes. The vast majority of the spacecraft's time will be spent under 
normal operation. As a result, the design of the spacecraft's Attitude Dynamics and 
Control Subsystem was driven by this mode. And, that is why a vast majority of the 
discussion will be on the Normal, On-Station operating mode. The other mode operations 
will largely feed from how the Normal/On-Station operations are conducted. [Ref. 6: p. 
340] 


6. Slew Mode 
We will utilize the Zero Momentum Three Axis Control System to simplify the 


solar array design. The spacecraft will continually slew about the unconstrained yaw axis 
in order to keep the solar array axis perpendicular to the solar radiation. Thus, the 
reduced array size possible with two degrees of freedom can be achieved with one array 
axis drive and continual spacecraft rotation about the yaw axis. It is important to ensure 
the reaction wheel selected can deliver the required torque. Slewing about the roll and 
pitch axis will not normally occur. If however, commanded to do so, the process and 
dynamics of slewing about the roll and pitch axis will be the same as that about the yaw 
axis. [Ref. 6: p. 351] 


7. Contingency / Safe Mode 

This mode is enabled in emergencies or when a certain level of uncertainty exists 
concerning the operational status of the spacecraft. In the Contingency/Safe Mode all 
non-essential equipment and actuators are turned off. Then, when given approval, the 
spacecraft will enter the Sun Acquisition Mode followed by extensive testing as the 


various components and actuators are brought back on line. 


8. Orbit Maintenance / Orbit Transition Mode 
If the spacecraft needs to be repositioned due to accumulated orbital error or a 


newly assigned orbit position it will enter the Orbit Maintenance/Orbit Transition Mode. 
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The first step in this mode is to shut down non-essential operations. Then, the spacecraft 
will orient itself such that the yaw axis of rotation points in the desired direction of 
movement. Then a spin up will take place followed by the firing of appropriate control 
thrusters. When the spacecraft has reached its new position satisfactorily, it will enter into 
the acquisition mode. 

Table 3-1 summarizes the type of spacecraft control and the sensors utilized by 


each attitude control mode. 
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MODE SENSORS ACTUATORS FUNCTION 


Orbit Insertion | - None - None Insert spacecraft 
into final orbit 


Acquisition - Sun Sensors (Pri) - Torque Rods (Pri) Determine 
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craft 
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- Star Trackers (Pri) 
- IRU (Sec) 










- Reorient 
Spacecraft as 
required 


Slew - Reaction Wheels (Pri) 
- Torque Rods (Sec) 


- Thrusters (Sec) 



































Contingency / | - Sun Sensor #1 (Pri) - Reaction Wheels (Pri) | Ensures spacecraft 
Safe - Earth/Sun Sensor #1 - Torque Rods (Sec) remains within 
(Pri) - Thrusters (Sec) limited operating 
- Sun Sensor #2 (Sec) constraints 
- Earth/Sun Sensor #2 







(Sec) 
- GPS (Pri) 
- IRU (Sec) 



















| Orbit 
Maintenance / 
Orbit 
Transition 
Mode 


spacecraft 
reposition / 
change in orbital 
parameters 


- Thrusters (Pri) 
- Thrusters (Sec) 





Table 3-1: Type of Spacecraft Control by Attitude Control Mode 


35 


C. THE MOTION OF THE SATELLITE AND THE SOLAR 
ARRAYS 


As a spacecraft decreases in altitude the time spent in eclipse becomes a 
considerable part of the total lifetime of the satellite. For instance, at 11,000 km the 
maximum eclipse will occur when the Sun is in the orbital plane and the eclipse time will 
be 12% of the orbital period. [Ref. 37: pp. 75-78] 

To take maximum advantage of the solar radiation, the solar arrays must be sun 
pointed. At the same time, the payload requires a constant nadir pointing configuration. 
At every position in the orbit these two demands must be satisfied. This will be 
accomplished with a combination of yaw steering and solar array pointing. Assuming the 
satellite has an arbitrary earth-pointing attitude, it must first rotate about the Z , -axis (axis 
parallel to the antenna; yaw-axis) with an angle wy, so that the axis about which the solar 
arrays rotate (BAPTA-axis) 1s in the plane perpendicular to the direction of the solar 
radiation. The second part of the motion consists of a rotation of the solar arrays about 
the BAPTA-axis with an angle n, so that the arrays point to the Sun. These two motions 
will now be looked at in much greater detail. [Ref. 37: pp. 75-78] 

The angles n and y are functions of time. They also depend on the position of the 
satellite in the orbital plane and the position of the Sun with respect to the orbital plane. 
To aid in our analysis the following definitions are made: First, 8B is defined as the angle in 
the orbital plane between the point with minimum distance (sub-solar point) to the Sun 
and the satellite, and 5 as the angle between the orbital plane and the direction of the Sun. 
This 1s.illustrated in Figure 3-2. Because of the direct relation between B and time, y and 


1 will be expressed as functions of 8 and 6. [Ref. 37: pp. 75-78] 
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LZ 


Figure 3-2: Definition of the angles yw, n, B and 6. 


The derivation of these functions y(B,5) and n(B,5) will be done below by means 
of vector algebra. [Ref. 37: pp. 75-78] 


The direction of satellite axes X;, Y >, Z, and of the sunlight s in the inertial 


XYZ coordinate system are: 


a7 


— sin(B)*sin(y) sin(B)*cos(y) —cos(f) 
Me =| cos(f)*sin(y) |; i =| cos( B)*cos(y) |; Ze =|-—sin(f) |; 
cos(y) — sin(y) 0 


—cos(d) 
sao) 
sin(d) 

While the vector of the sunlight s stands perpendicular to the BAPTA-axis, and 
because the BAPTA-axis and the Y .-axis are parallel, Y. is also perpendicular to s. 
Now w(B,5) can be determined: 
¥. -5 =—sin(P) cos(w) cos(d) — sin(y) sin(d) = 0 

~sin(f) ). — sin(f) 


a tn(y)=( tan(d) } ae am un“ eno 


The direction of the solar panels inthe X., Y., Z. satellite coordinate system 


— cos(77) 
(see figure III.2) is: ps = 0 


sin(77) 


The direction of the solar panels in the XYZ system is obtained by performing the 


matrix multiplication p = Bs nove Zs |*[Bs|: 


sin(B)*sin(y )*cos(7) — cos(8)*sin(7) 
=> p=i-—cos(f)*sin(y)*cos(7) — sin(Z)*sin(77) 
— cos(y )* cos(77) 


Since p and s have to point exactly in opposite direction, also 7(,0) can be 


determined: 
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x y 2 
pxs=det|p, Pp, p;}= 
5, 5 3S; 


(— cos(f)* sin(y)* cos(77) — sin( 8)* sin( 7))* sin(d) 
(- sin( f)*sin(w)*cos(77) — cos(8)*sin( ™))* sin(d) + cos(w)*cos(7)*cos(d) | = 0 
— (cos( f£)*sin(y)*cos(77) + sin( B)*sin( 7))* cos(d) 





__siny) 8, ae) 
an 


w(B) and n(f) are plotted in Figures 3-3 and 3-4. From these figures it is clear 
that when 6 — 0° (when the Sun is almost in the orbital plane) the satellite has to make 
radical short-time movements. Especially the motion of the satellite about the yaw axis 
nearby B = 0° and B = 180° is critical because of the relatively large moment of inertia 
about this axis (I _ ) and the high rotation speed (see Appendix C). [Ref. 37: pp. 75-78] 

For this reason a limit value 6, was chosen. If the satellite is in an orbit with 6 > 
8 | then an alternative motion will be executed. The yaw-axis will keep a constant 
attitude: w = 90°, so that the BAPTA-axis is perpendicular to the orbital plane. In this 
case the solar arrays would rotate 360°/orbit in one direction. The power lnes would 
therefore turn around the BAPTA-axis and after a while make the rotation of the arrays 
impossible. The only way to anticipate this is rotating them in the opposite direction 


during the eclipse. This alternative motion is plotted in Figure 3-5. [Ref. 37: pp. 75-78] 
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Figure 3-3: Motion of Yaw-axis (y) and solar array (yn) during one orbit, starting at the 


sub-solar point (1.e., Noon) in terms of the sun rotation angle. 
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Figure 3-4: Motion of Yaw-axis (yw) and solar arrays (n) during one orbit, starting at the 


sub-solar point in terms of the time in minutes. 
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Figure 3-5: Alternative motion of yaw-axis (yw) and solar arrays (1) during one orbit, 


starting at the sub-solar point, for |6| > 5°. 
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IV. QUANTIFY DISTURBANCE ENVIRONMENT 


In this part of the design we must understand what disturbances are acting on the 
spacecraft and then try to quantify them. Disturbance torques can be either external or 
internal. External disturbance torques consist of gravity gradient effects, magnetic field 
torques acting on the vehicle, the effects of solar radiation, and aerodynamic effects. 
These external disturbance torques are most likely the largest source of disturbance to the 
spacecraft and may be easier to quantify. Internal disturbance torques come from within 
the spacecraft. As a result, we have a great amount of control over these. [Ref. 6: p. 352] 

In this section an investigation of both external and internal disturbance torques 
will be performed. And, a knowledge of both magnitude and direction of disturbance 
torques and angular momentum will be attamed. This information will be very useful in 


selecting and sizing the ADCS hardware. [Ref. 6: p. 352] 


A. EXTERNAL DISTURBANCE TORQUES 


From the standpoint of stability and control, the principal external torques come 
from four environmental effects: solar radiation pressure, gravity gradients, aerodynamics 
and magnetics. Solar radiation pressure is generally a significant source of attitude and 
trajectory errors for high altitude (> 1000 km). Gravity gradients, which result from the 
extended dimensions of the spacecraft, may either cause disturbing or perturbing torques 
or provide restoring torques when the effect is used for attitude control. Gravity 
gradients, as well as the magnetic torques caused by the interaction of the spacecraft's 
magnetic materials with the planetary magnetic field are most significant at low altitudes 
(< 1000 km). Similarly, aerodynamic effects are significant only below 500 km altitude 
and are generally negligible above 1000 km. Finally, aerodynamic torques acting on the 
spacecraft are functions of spacecraft geometry, attitude and altitude. Figure 4-1 


illustrates the change in external disturbance as a function of altitude. [Ref. 22: pp. 77-78] 
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Figure 4-1: Typical torques on a satellite as a function of altitude 


1. Solar Pressure Torques 


The solar pressure torque is the major long-term disturbance torque for a 
spacecraft at higher altitudes. The solar radiation forces are due to photons impinging on 


the spacecraft surfaces. In general a fraction, p,, of the impinging photons will be 
specularly reflected, a fraction, p,, will be diffusely reflected, and a fraction, p,, will be 
absorbed by the surface (see Figure 4-2). The solar pressure moment M, is given by 


equation 4-1. [Ref. 25: pp. 133-134] 


ge 
Surface Area = A Specularly 


Solar Pressure = P 





reflected 
photons <a. as 


- 
—_ 


—_ 


Incoming <——————————_ 


$ 
photons 
a <—_—_____—_ 
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Figure 4-2: Solar Radiation Force on a Spacecraft Surface 


M. = PA(n-S)F x (1-2) a p, to) (4-1) 


where F is the vector from the center of mass of the spacecraft to the center of pressure of 


a given area, A. The solar radiation pressure, P , is generally assumed to be constant and 


to have the value 4-644x1° N/m’. [Ref. 25: pp. 133-134] 


2) Gravity Gradient Torques 

A spacecraft body experiences a gravity-gradient torque due to the variation of 
distances between the spacecraft mass points and the center of the mass of the earth. The 
gravity gradient has been used on early low earth-orbit satellites to maintain the earth 
pomting of antennas or other instruments. The gravitational torques about the principal 


body axes e;, e2, and e3 are given in equation 4-2. [Ref. 25: p. 131] 
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T;? = —K(1,, — 1,)cos(@,)sin(@;)sin(O;)) =" fons Gao. 0 
T}® = K(/,, — 1,3) cos(6;) cos(@,) sin(9, ) (4-2) 
T{® =~K(1,,- I,,)cos(6;)[cos(0, yi sin(8, ). 


3. Aerodynamic Torques 


The satellite will, in general, pass through an atmosphere of density p, with a 


velocity ¥. The magnitude of the aerodynamic force F‘” is then given as 
(nee 
F =~ pi- FAC, (4-3) 


where A is the reference area of the satellite (such as the cross section along v) and C, 


is the total drag coefficient. Torque is given by 


l 


io i veiSG (4-4) 
or 
T°) = = pC,Av'(ax/,) (4-5) 


where / is the length of the perpendicular from the mass center to the force line of action 
and w is the unit vector in the velocity direction. [Ref. 22: p. 90] 

Difficulties arise in determining p,/,and C,. At satellite altitudes, o is highly 
dependent on the time of day and the level of solar activity. For example, at 600 km the 
solar daytime maximum density and the nighttime minimum density may differ by a factor 
of 100. Table 4-1 gives typical values of the daytime maximum air density as a function of 
orbital height. It is important to note that as altitude increase the aerodynamic pressure 
effects decrease. As a result, our satellite orbit of 11,000 km will have negligible 


aerodynamic pressure effects. [Ref. 22: p. 90] 
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Table 4-1: Density as a function of height above the earth's surface 





For a satellite having a spherical shape, an average value of C, = 2.2 can be taken 
which is computed assuming "free molecular flow" (i.e., the molecular mean free path is 
assumed to be large compared to the size of the satellite and, therefore, inter-particle 


collisions are ignored). For a cylinder C, may be taken as 3. [Ref. 22: p. 90] 
4. Disturbance Torques Due to the Geomagnetic Field 


a. Geomagnetic Field 

The geomagnetic field can be represented to a first order approximation as 
a magnetic dipole with the axis inclined to the Earth's spin axis by approximately 11.5°. 
Observations have shown that the sources of the geomagnetic field are in the core and at 
the surface of the Earth as well as in the upper atmosphere. The primary source is 
presumed to be a system of electric currents in the Earth's molten interior. The changes in 
these currents account for the migration of the geomagnetic poles on the Earth's surface. 
[Ref. 22: p. 84] 

The geomagnetic field, B, is a function of altitude, longitude, and latitude. 
And, can be expressed in units of gauss, gamma (y), or tesla (webers per square meter). 
A one gauss magnetic field intensity acting on a unit magnetic pole produces a force on 


one dyne. The magnitude of the field is given by equation 4-6. [Ref. 22: p. 84] 
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B 








= =~ +3sin oy (4-6) 


where @,, is the magnetic latitude measured from the geomagnetic equator. At the earth's 
surface along the geomagnetic equator, 0,,= 0, and B, = 0.311 gauss. At the magnetic 
pole, B, = 0.622 gauss or twice the equatorial value. [Ref. 22: p. 84] 


A plot of the Earth's equatorial magnetic field intensity B,, in tesla units as 


a function of altitude is shown in Figure 4-3. The effect of latitude is given in Figure 4-4, 


where for any latitude B, = @B,,. [Ref. 22: pp. 87-88] 
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Figure 4-3: Earth's Magnetic field intensity vs altitude above the magnetic equator. 
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Figure 4-4: Multiplication factor a for magnetic field intensity as a function of magnetic 


latitude. For any latitude B, = 


equatorial ° 


b. Magnetic Torques 
Magnetic torques acting on a spacecraft can result from the interaction of 
the spacecraft's residual magnetic field and the geomagnetic field. Thus, if M is the sum 


of all magnetic moments in the spacecraft the torque acting on the spacecraft is 


T” =MxB (4-7) 


where B = the geomagnetic field vector. In general, M can be caused by permanent and 
induced magnetism or by spacecraft-generated current loops. This is how Magnetic 
torque rods operate and can be used for momentum dumping. The units of / may be 
gauss-cm*, ampere-m? , or pole-cm. For example, if M is in ampere-m’, then T°” is in 


Newton-meters. [Ref. 22: p. 89] 
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The concept of Magnetic Torquers is simple, but there is a difficulty in 
quantifying external disturbance torques due to the geomagnetic field. M is the magnetic 
field produced by the spacecraft and 1s generally unknown until the spacecraft is built and 
tested. Once tested, an estimate for the magnetic disturbance torque can be made. [Ref. 


22 bo 


B. MATLAB SIMULATION OF EXTERNAL DISTURBANCE 
TORQUES 


To attain an accurate quantifiable value for external disturbance torques a matlab 
simulation was developed for each of the external disturbances. A brief explanation of 


each scenario follows: 


1. Solar Pressure Torques 

1. Developed the Simulink model appropriate for the Solar Pressure Simulation 
(see Appendix E). 

2. Established the initial conditions for the velocity and position vector of the 
spacecraft relative to the earth centered inertial system (see Appendix E). 

3. Defined various constants and solved for v and 6 as a function of time. v is the 
earth centered true anomaly and 6 is the sun orbital element analogous to v (see Appendix 
Ey 

4. Established numerous direction cosine matrices. This part of the process was 
very tedious. The details are excluded because it is assumed the reader is familiar with the 
various orbit frames of reference (see Appendix E). 

5. Computed the unit vector from the spacecraft to the sun (see Appendix E). 

6. Defined unit vectors that are perpendicular to each exterior piece of the 
spacecraft (see Appendix E). 

7. Defined dr vectors from the center of mass to each external piece's center of 
pressure (see Appendix E). 


8. Performed various vector operations (see Appendix E). 
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9. Calculated the moment due to solar pressure as given by equation 4-1 [Ref. 25: 
p. 133]: 


M, = PA(neS)F x ( + p,)S+2(p,+ (2) ba 


where, 


M,= The solar radiation pressure moment. 


r= The vector from the center of mass of the spacecraft to the 
center of pressure of a given area, A. 


P= The solar radiation pressure. 
A= Area 


Le 


Fraction of photons impinging on the spacecraft's surface 
that are specularly reflected. 


P, = Fraction of photons impinging on the spacecraft's surface 
that are diffusely reflected. 


10. Computed the M,, M,, and M, for solar pressure torques for each instant in 
time (see Appendix A-10 to A-17). 


11. Solved for angular velocities iteratively (see Appendix E). 


12. Utilized Poisson's equations to solve for the dem components from the body 
frame to the orbit frame iteratively (see Appendix E). 


13. Solved for rand v iteratively utilizing state space variables (see Appendix E). 


Figures 4-5 through 4-7 are the result of the calculations described above. A brief 
explanation of each figure follows: 
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a. Figure 4-5 


Figure 4-5 is a graph of the three components of torque, (i.e.. 


T,, T,, and T,) acting on the spacecraft due only to solar pressure disturbances, when the 


spacecraft is not being controlled by the Attitude Dynamics and Control Subsystem. Also, 


the torque components are with respect to the spacecraft's body fixed frame. 


b. Figure 4-6 


Figure 4-6 is a graph of the three components of the spacecraft's angular 
velocity due only to solar pressure disturbances when the spacecraft is not being 
controlled by the Attitude Dynamics and Control Subsystem. Also, the angular velocity 


components are expressed with respect to the spacecraft's body fixed frame. 


c. Figure 4-7 


Figure 4-7 is a graph of the three components of torque, (i.e., 


T., T,, and T,) acting on the spacecraft due only to solar pressure disturbances, when the 


spacecraft is being controlled by the Attitude Dynamics and Control Subsystem and is 
maintaining a nadir pointing configuration. Also, the torque components are with respect 
to the spacecraft body fixed frame. 

It should be noted that the solar radiation torque component along the 
pitch axis is periodic and secular, with a period equal to the orbit period. Hence the net 
effect of this torque is a build up of Angular Momentum. On the other hand, the 
components along the other axes have cyclic components, resulting in no net change in the 


spacecraft angular momentum over a full orbit. 
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Figure 4-5: The three components of torque due to solar pressure as a function of time 
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with no control. 
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Figure 4-6: The three body components of angular velocity due to solar pressure as a 


function of time, with no control. 
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Figure 4-7: The three components of torque due to solar pressure as a function of 


time, with active control. 


2. Gravity Gradient Torques 

1. Developed the simulink diagram (see Appendix E). 

2. Defined constants and user friendly conversions (see Appendix E). 

3. Determined initial conditions for r, v , and the dcm between orbit frame and 
principal inertia frame (see Appendix E). 

4. Solved for angular velocities iteratively (see Appendix E). 

5. Utilized Poisson's equations to solve for the dcm components from the body 
frame to the orbit frame iteratively (see Appendix E). 

6. Solved for r and v iteratively utilizing state space variables (see Appendix E). 


7. Calculated the moment due to Gravity Gradient as given by equation 4-11. 


de = ie = La 
Le = K(1,, a Eanes 
2 = Ky ede aga 
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Figures 4-8 and 4-9 are the result of the calculations described above. A brief 


explanation of each figure follows: 


a. Figure 4-8 


Figure 4-8 is a graph of the three components of torque, (i.e. 


I, 1,, and T,) acting on the spacecraft due only to gravity gradient disturbances, when 


the spacecraft is not being controlled by the Attitude Dynamics and Control Subsystem. 
Also, the torque components are with respect to the spacecraft's body fixed frame and 


increase over time such that they cause the spacecraft to become unstable. 


b. Figure 4-9 


Figure 4-9 is a graph of the three components of the spacecraft's angular 
velocity due only to gravity gradient disturbances when the spacecraft is not being 
controlled by the Attitude Dynamics and Control Subsystem. Also, the angular velocity 


components are expressed with respect to the spacecraft's body fixed frame. 
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Figure 4-8: The three components of torque due to gravity gradients as a function of 


time, with no control. 
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Figure 4-9: The three body components of angular velocity due to gravity gradients as a 


function of time, with no control. 


3. Aerodynamic Pressure Torques 


1. Calculating the Aerodynamic Pressure Torques was very similar to calculating 
solar pressure torques. Nearly the same technique was used except for the following (see 
Appendix E). 


2. The fundamental equation used to compute is given by equation 4-5 [Ref. 6: p. 
309]: 


M = (4) pc,av (i, x lias} 


where, M = Torque due to Aerodynamic Pressure 
p = Atmospheric density 
C, = Coefficient of drag 


A= Area perpendicular to w, 


u, = Unit vector in velocity direction 


ie = Vector distance from center of mass to center-of- 
pressure 

3. Computed the density (p) as a function of height above the earth (see 
Appendix E). 

4. Computed the cross sectional areas of the spacecraft in the three body frame 
coordinate axis directions (see Appendix E). 

5. Used the same Euler rates and Poisson's equations to solve iteratively (see 
Appendix E). 

6. Utilized state space variables to solve for rand v iteratively (see Appendix E). 

Figures 4-10 through 4-12 are the result of the calculations described above. A 


brief explanation of each figure follows: 


a. Figure 4-10 


Figure 4-10 is a graph of the three components of torque, (i.e., 


T., 1,, and T,) acting on the spacecraft due only to aerodynamic pressure disturbances, 


when the spacecraft is not being controlled by the Attitude Dynamics and Control 
Subsystem. The torque components are with respect to the spacecraft's body fixed frame. 
Also, these torque magnitudes are orders of magnitude less than gravity gradient and solar 
pressure disturbances. And, catastrophic cancellation (i.e., computer noise) most likely 


occurred when calculating 7,. This was probably due to the fact that the torque values 


were so small. Asa result, the torques due to aerodynamic pressure are inconsequential. 


b. Figure 4-11 

Figure 4-11 is a graph of the three components of the spacecraft's angular 
velocity due only to aerodynamic pressure disturbances when the spacecraft is not being 
controlled by the Attitude Dynamics and Control Subsystem. The angular velocity 


components are expressed with respect to the spacecraft's body fixed frame. Also, the 


57 


aerodynamic torques are orders of magnitude less than gravity gradient and solar pressure 


disturbances. Therefore, the effects are inconsequential. 


c. Figure 4-12 


Figures 4-12 is a graph of the three components of torque, (i.e., 


I, [,, and T,) acting on the spacecraft due only to aerodynamic pressure disturbances, 


when the spacecraft is being controlled by the Attitude Dynamics and Control Subsystem 
and maintaining a nadir pointing configuration. The torque components are with respect 
to the spacecraft body fixed frame. Also, the effects of aerodynamic pressure disturbances 
are orders of magnitude less than gravity gradient and solar pressure disturbances. Also, 
catastrophic cancellation (i.e., computer noise) may have occurred. As a result, the 


torques due to aerodynamic pressure are inconsequential. 
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Figure 4-10: The three components of torque due to Aerodynamic Pressure as a function 


of time, with no control. 
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Figure 4-11: The three body components of angular velocity due to Aerodynamic 


Pressure as a function of time, with no control. 
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Figure 4-12: The three components of torque due to aerodynamic pressure 


disturbances as a function of time, with active control. 


4. Combined External Disturbance Torques 

Next, calculations were made for when all three torques acted on the spacecraft 
simultaneously. This required combining the three separate Matlab programs into one 
program. The motion of the spacecraft was then simulated and the sum of the 
environmental disturbance torques computed (see Appendix E). 

Finally, an estimation of the total torque imparted to the spacecraft over one year 
for each of the four situations was performed. To do this for gravity gradient and 
aerodynamic required taking data from several orbits, then computing the average torque 
over that period and finally multiplying by the time expired. This method is good for 
gravity gradient and aerodynamic as long as data from several orbits is computed. 
Computing the total torque from the solar pressure is more complicated. Since the Beta 
angle changes with time throughout the year, several orbits will not be adequate. So, the 
simulation was ran for a delta t of one year, then the total torque was computed. The 


results can be seen in Table 4-2. 


SOLAR GRAVITY AERO- COMBINED 
PRESSURE GRADIENT | DYNAMIC 
Average 7, 9.33e-10 Nm | 20 | 20 | -2.5e-6 Nm 


Average T, 1.33e-6 Nm 1.2e-6 Nm 9.55e-29 Nm ].4e-4 Nm 
Average T, 5.03e-8 Nm 6.6e-11 Nm | -2.78e-6 Nm 











Maximum Torque 
C2’) 1.8e-6 Nm le-4 Nm Je-11 Nm 2.2e-4 Nm 


0 
0 , 
Estimated H, (1 year) | 0.0296 Nms 0.1718 Nms 
Estimated H, (1 year 42.23 Nms 3155 Nms_ | 3.03e-21 Nms | 3197.23 Nms 
1.594 Nms 0.0021 Nms | 0.19033 Nms 


Table 4-2: Summary of External Disturbance Torques and Angular Momentum Imparted 
to the Spacecraft 


Estimated H, (1 year 


Figures 4-13 through 4-15 are the result of the calculations described above. A 


brief explanation of each figure follows: 
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a. Figure 4-13 
Figure 4-13 is a graph of the three components of torque, (i.e. 


T,, T,, and T,) acting on the spacecraft due to solar pressure, gravity gradient and 


aerodynamic disturbances, when the spacecraft is not being controlled by the Attitude 
Dynamics and Control Subsystem. Also, the torque components are with respect to the 
spacecraft's body fixed frame. 


b. Figure 4-14 


Figure 4-14 is a graph of the three components of the spacecraft's angular 
velocity due to solar pressure, gravity gradient and aerodynamic pressure disturbances 
when the spacecraft is not being controlled by the Attitude Dynamics and Control 
Subsystem. Also, the angular velocity components are expressed with respect to the 


spacecraft's body fixed frame. 


c. Figure 4-15 


Figures 4-15 is a graph of the three components of torque, (i.e. 


T,, T,, and T,) acting on the spacecraft due to solar pressure, gravity gradient and 


aerodynamic pressure disturbances, when the spacecraft is being controlled by the Attitude 
Dynamics and Control Subsystem. Also, the torque components are with respect to the 
spacecraft body fixed frame. 

It should be noted that the total torque component along the pitch axis is 
periodic and secular, with a period equal to the orbit period. Hence the net effect of this 
torque is a build up of Angular Momentum. On the other hand, the components along the 
other axes have cyclic components, resulting in no net change in the spacecraft angular 


momentum over a full orbit. 
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Figure 4-13: The three components of torque due to total external disturbance torques as 


a function of time, with no control. 
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Figure 4-14: The three body components of angular velocity due to total external 


disturbance torques as a function of time, with no control. 
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Figure 4-15: The three components of torque due to solar pressure, gravity 


gradients and aerodynamic pressure disturbances as a function of time, with active control. 


C. INTERNAL DISTURBANCE TORQUES 

There are several other kinds of disturbance torques that are internal to the 
spacecraft. First, these torques may be due to moving parts such as rotating wheels, 
circulating fluids or tracking devices. Fortunately, the ADCS designer has some control 
over these torques. For example, if a disturbance torque in one component is found to be 
much larger than the others, the component can be changed in order to meet higher 
specifications. However, there is a trade-off here. A change can be made to meet higher 
standards only at the expense of more mass and/or cost. Next, an internal torque may 
arise due to thruster mis-alignment with the center of gravity. This will only occur during 
thruster firing and will be corrected in a closed loop Attitude Dynamics and Control 
Subsystem. The concern here is that the disturbance will not cause the spacecraft to 
exceed the minimum pointing requirements. [Ref. 6: p. 352] 

Of greater importance is slosh and mechanical machinery torques. They can, 
however, be controlled by proper selection of hardware. Table 4-3 summarizes the 


common internal disturbance torques. [Ref. 6: p.352] 
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Effect on Vehicle Typical Values 


Uncertainty in center of e Unbalanced torques during | | to 3 cm 
gravity (cg) firing of coupled thrusters 

e Unwanted torques during 
translation thrusting 


Thruster misalignment Onto Urs 


Mismatch of thruster Similar to cg uncertainty + 5% 
outputs 


Rotating machinery (pumps, | Torques that perturb both Dependent on spacecraft 
tape recorders) stability & accuracy design; may be 

compensated by counter- 
rotating elements 

































































Liquid sloshing Torques due to fluid motion Dependent on specific 
and variation in center of mass | design; may be controlled 
location by bladders or baffles 

Dynamics of flexible bodies | Oscillatory resonance at Depends on spacecraft 






bending frequencies, limiting structure 
control bandwidth 
Attitude disturbances when 
entering/leaving eclipse 

























Thermal shocks on flexible 
appendages 


Depends on spacecraft 
structure. Worst for 

gravity gradient systems 
with long inertia booms 
Table 4-3: Principal Internal Disturbance Torques which spacecraft designers have a great 





deal of control over. 
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V. ATTITUDE DYNAMICS AND CONTROL SUBSYSTEM 
HARDWARE SELECTION AND DESCRIPTION 


A. SIZING THE ADCS REACTION WHEELS 


1. Torque and Angular Momentum Requirements Due to Yaw 
Steering 
In Chapter III the yaw steering requirements of the spacecraft were investigated 


and it was concluded the Yaw command angle could be stated as: 





= _,| —sin(Z) : 
y(f,0) = tan tan(S) (5-1) 


where, 


B is the angle in the orbital plane between the point with minimum distance (sub- 
solar point) to the Sun and the satellite. And, 6 is the angle between the orbital plane and 
the direction of the Sun (see Figure 5-1). 

By looking at Figure 3-4 it is apparent that the angular velocity and angular 
acceleration maximums occur for the minimum value of 6. And, since our minimum 6 is 
equal to 5° we will single this curve out and find the maximum angular velocity and 
maximum angular acceleration the spacecraft must perform in the yaw direction. Since 
this requirement is the most demanding; if we can show it will be met then any other 


torque and angular momentum requirement will also be satisfied. 





Figure 5-1 is a graphof w(,5°) = an] for 6 = 5°. It can readily be 


shown that y (f,5°) which is the angular velocity is given by 


Te 01 506) | ens 6e= 4 =) | - 


+ 130.646 * [sin(2.76e — 4* 1) 


and, yw (f,5°) which is the angular acceleration is given by 
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sin(2.76e —4* 1) | 


CCS): = 
ck 1 + 130.646 *|sin(2.76e — 4 *1)] 


(5-3) 


Por [cos(2.76¢ — 4 *1)] * sin(2.756e — 4* 1) 


1 + (130.646) *[sin(2.76¢ —4 * of] 


Figures 5-2 and 5-3 are graphs of the angular velocity and angular acceleration 


respectively for 6 = 5°. 





Figure 5-1: Graph of yw(Z,5°), where y is in degrees and time is in seconds. 
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Figure 5-2: Angular Velocity in units of degrees/second. 
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Figure 5-3: Angular Acceleration in units of degrees/second° . 
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From the graphs it is observed that the maximum angular velocity occurs at t=0, 
11400 and 22800 seconds and the maximum angular acceleration occurs at t = 0°, 
11400° , 11400° , and 22800 seconds. The respective values are as follows: 


y (B85). = 0.18 / sec = 3.142¢ — 3 rad / sec 
wy (B.S) n» * 0.00037 / sec” = 6.4577e — 6 rad / sec’ 


The maximum acceleration of the satellite has to be realized by accelerating the 
reaction wheel in the opposite direction. The maximum torque T _,, needed follows from 


the formula: 


To. = Lop*W max = (14113) *(6.4577e — 6 rad / sec”) = 0.0911 Nm (5-4) 


m 


And, the maximum angular momentum B __.. is given by: 


B= Ioz*W max = (14113)*(3.142e — 3 rad / sec) = 44.34 Nms (5-5) 


The total acceleration Aw of the reaction wheel due to the rotation can simply be 


calculated with the following relation: 


= wt (5-6) 


@ 


Ly*W mx =1,*A@ => A@ 


where, we find J, from the formula for the moment of inertia for a hollow cylinder 
with a relative small wall thickness, J = mR*. So, for a reaction wheel with a mass of 


10.5 kg and a radius of 0.1965m: 


=> 1,~04054 kgm’ 


such that, 


[ a 
A@ = ee = 9.69 rad / sec = 93 rpm 


@ 
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In conclusion, the reaction wheel to be selected must meet these requirements for 


torque, momentum and acceleration. 


2. Torque and Angular Momentum Requirements Due to External 
Disturbances 
The maximum torque and angular momentum requirements due to external 


disturbances were calculated and found to be: 


Ti max = 8€—-4 Nm FA, na, = 9.1718 Ns 
D max = 2-2€—-3 Nm FA, ax = 9-498 Nms 
Ty max = OC —4 Nm FT vax = 9.19033 Nms 


As can be seen, these torque and angular momentum requirements are not as 
demanding as those resulting from yaw steering. However, the external disturbances 
produce a secular torque which builds up angular momentum with time. As can be seen 
from Figure 4-16; a secular torque results in the 7, direction and creates a buildup of 


angular momentum as given by 


(7; 


yng ) X (1 month) = 362 Nms 

over the course of one month. This build up must be dumped by either the Magnetic 
Acquisition / Despin system or the spacecraft thrusters. Not accounting for these secular 
disturbances appropriately could result in a premature mission due to unexpected demands 


of spacecraft fuel. 


B. SIZING THE MAGNETIC TORQUERS 
Magnetic torquers use the Earth's magnetic field, B , and electrical current through 
the torquer to create a magnetic dipole ( D) that results in torque (7 ) on the vehicle given 
by: 
T 
ee 5-7 
F (5-7) 


Magnetic torquers used for momentum dumping must equal the peak disturbance 
plus a margin to compensate for the lack of complete directional control. 
From equation 5-8 we can estimate the worst case earth magnetic field for our 


orbit given by using the following values: 
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ne |i aeeee me 
go 
Uy, =8.1e25 g-cm’ =8lelST-m 


B 


Oo 








= ££(143sin0,)"” (5-8) 
r 
= IBnax{ = 155e-6 T 


Next, we calculate the torque rod's magnetic torquing ability (dipole) to counteract 
the worst case external disturbance; 


T ax = 2-2e-3 Nm 


ie 


D=—=1419.35 A-m° 
B 


which is achievable. The earth's field is cyclic at twice the orbital frequency; thus, 
maximum torque is available only twice per orbit. A torquer of 2500 to 3000 
A-m capacity will provide sufficient margin to counteract these expected maximum 
disturbance torques and/or provide a momentum dumping capacity for the reaction 
wheels. 


C. OTHER SENSOR SELECTION ISSUES 

The driving requirement for the remaining sensors and actuators are the attitude 
knowledge and pointing requirements as stated in Chapter 6 under the system level 
requirements section. All sensors and actuators selected meet these requirements with a 


considerable margin built in. 
D. HARDWARE SELECTION 


1. Spacecraft Control Processor 

The Satellite Control Processor unit selected is produced by Honeywell and is 
scheduled to be used on the (International Space Station. Honeywell's 
Multiplexer/Demultiplexer (MDM) is a high speed interface and computer that controls 
sensors, effectors and other unique devices. [Ref. 29] 
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The MDM collects signals from sensors, effectors and subsystems; then processes 
the information and transmits signals and commands. For example, the MDM directs the 
solar arrays to point at the sun. [Ref. 29] 

The MDM 1s rugged and adaptable to a wide range of applications. Five standard 
I/O circuit cards cover most applications, but this can be modified if needed. Also, the 
MDM size can be modified based on the number of I/O circuit cards needed. [Ref. 29] 

To enhance reliability and supportability, software and hardware self-test features 
can detect and report 94% of all possible failure modes. In addition Honeywell can help 
develop I/O cards and user application programs that meet individual needs. Figure V-4 is 
an overview of the MDM architecture. Features of the MDM: [Ref. 29] 


1. The basic MDM consists of a chassis, power supply, I/O Control Unit, I/O 
backplane, SX backplane and I/O circuit cards selected to meet specific requirements. 

2. The I/O Control Unit includes the processor, memories, 1553 bus interface 
function, I/O Circuit Card Control logic, and Analog-to-Digital converter. 

3. Communication is via independent, dual redundant MIL-STD-1553B channels, 
enabling the MDM to operate in Bus Controller, Bus Monitor or Remote Terminal mode. 
Firmware supports functions of each 1553 channel independently of the user application 
software. The MDM is capable of serial asynchronous communication via RS449 with 
RS485 driver/receivers and provides 16 bits of parallel I/O functionality. The MDM self- 
test capabilities include Built-In Test (BIT) software routines and supporting hardware 
Built-In Test Equipment (BITE). 
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Figure 5-4: MDM Architecture 
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2. Magnetic Acquisition and De-spin System 

The Magnetic Acquisition and Despin (MAD) system is designed for use on 
spacecraft for the reduction of spin rates and momentum dumping and is built and 
manufactured by ITHACO Space Systems. The spin rates can be encountered as a 
normal part of a mission or can result from ground control errors or onboard equipment 
failures. Under normal conditions, the launch vehicle leaves the satellite spinning at the 
end of orbit insertion. If the spacecraft is to be stabilized, the spin rate imparted by the 
launch vehicle must be corrected. The unexpected activation of a thruster or a momentum 
wheel failure would also result in unacceptable spin rates. Corrective action can often be 
taken if the vehicle rates can be reduced in some manner. The Magnetic 
Acquisition/Despin system performs this function. [Ref. 30, Ref. 31] 

The MAD system consists of three Torque rods, an electronics controller, and a 
three-axis magnetometer. One Torque rod is aligned with each axis of the spacecraft and 
is mounted away from any instruments sensitive to magnetic fields. In addition, Torque 
rods are separated from one another in order to avoid cross coupling effects. The use of 
Torque rods allows complete momentum management of the spacecraft without 
consumables in almost any orbit, thus providing a clean, contamination free environment. 
Even in high orbits, Torque rods can make a significant contribution to the total 
momentum management problem. By eliminating or reducing the need for consumables, 
the use of Torque rods can reduce mission costs significantly. Next, the electronics 
controller can be mounted at any convenient location. The electronics should be 
accessible, however, since a trimming operation at the spacecraft level is required to 
compensate for the unavoidable magnetometer to torque rod coupling. The 
magnetometer should be mounted so that one Torque rod is aligned with each of the 
magnetometer axes. The magnetometer probe should be mounted as far away from the 
Torque rods as possible to minimize magnetic coupling. The maximum rate from which 
the MAD system can achieve capture depends upon the bandwidth that can be used. The 
maximum bandwidth in any given system depends on the amount of coupling between the 
Torque rods and the magnetometer. Thus, maximum separation is an advantage for the 


available capture capability of the MAD system. [Ref. 30, Ref. 31] 
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The MAD system relies upon a "minus B-dot" control law. Each Torque rod is 
driven to a value determined by the negative rate of change of the magnetic field along the 
rod. The effect is that the Torque rods behave like magnetic brakes and lock onto the 
Earth's magnetic grid. From this terminal state, the spacecraft's attitude control system 
can capture control or excess momentum can be dumped. The system specifications are 


as follows: 


pose S| wei | PowER 
0.9 ks N/A 


Table 5-1: Specifications of Magnetic Acquisition / Despin System 










3. Earth Sensor / Sun Sensor Assembly 


a. Overview 

The Earth Sensor/Sun Sensor Assembly (ES/SSA) is an innovative 
combined Earth horizon and two axis sun sensor which offers unprecedented performance 
in terms of accuracy, mass, power consumption, reliability, and automation. This sensor 
autonomously determines three-axis attitude information at medium earth orbits, 
geosynchronous, or super geosynchronous altitudes, as well as during transfer orbit. [Ref. 
32| 


b. Theory of Operation 

The Earth Sensor portion of the ES/SSA is a scanning sensor with two sets 
of two infrared (IR) Fields of View which scan about the axis of rotation (scan axis) at a 
rate of 2 revolutions per second. Each of the four IR scans possesses a 1.5° instantaneous 
Field of View. As the scan progresses, each of the instantaneous Fields of View crosses 
the relative "cold" of space to the "warmth" of Earth and back into space. Roll attitude is 
derived from differences in the sensed Earth chord widths. Pitch attitude is derived from 
the phase angle of the sensed Earth crossings. [Ref. 32] 

The sensor has also been designed to operate in transfer orbit and/or 
constant altitudes ranging from 11,000 km through 77,000 km. At geosynchronous 
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altitudes all four Earth scans cross the Earth's disc at nominally 3 and 6 degrees above the 
equator and 3 and 6 degrees below the equator. [Ref. 32] 

The Sun Sensor portion of the ES/SSA incorporates two static, mutually 
orthogonal sensing (axes) arrays to detect the orientation of the sun. Each axis 1s 
comprised of a slit allowing a plane of light to pass through and intersect with an 
orthogonally mounted, high resolution, linear photodiode array. The point of intersection, 
which indicates sun angle, is obtained by reading the energy stored in the array. [Ref. 32] 

Each ES/SSA is comprised of a scanning Earth sensor and a two-channel 
sun sensor housed in a single aluminum enclosure. Mounting is achieved by an 
asymmetric three point mount which allows for easy alignment and shimming if required. 
[Ref. 32] 

The Earth sensor operates over an attitude range of + 24° in pitch and + 
14° in roll, while the Sun Sensor typically operates with sun locations of + 60° in azimuth 
and + 30° in elevation, although fields as large as 128° x 128° may be accommodated. All 
processing and calibration data utilized to derive pitch and roll output from the Earth 
Sensor, as well as sun azimuth and elevation from the Sun Sensor is available on-board the 
sensor. Output occurs over a digital interface governed by the MIL-STD-1553 data bus, 
EIA RS-422A, or user specified interface. Upon asserting a data request the ES/SSA 
outputs all Earth Sensor and Sun Sensor information including diagnostics and status 
information. Table 5-2 lists the Earth Sensor / Sun Sensor Assembly Specifications. [Ref. 
32] 
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EARTH SENSOR 
Roll/Pitch Accuracy (Includes all long 
term & diurnally varying errors) + 0.02° (30) 
Transfer Orbit Roll / Pitch Accurac + 0.1° (30) 
Attitude Range + 14° Roll, + 24° Pitch 
SUN SENSOR 
Range Upto TZ23eal2 se 


Accuracy + 0.01° (30) 
OVERALL INSTRUMENT 


Manufacturer ITHACO Space Systems 
Power Consumption 
Mass 5.9) IDG) Se) 


Dimensions 6.4 x 3.9 x 8.0 in 
Output Data Rate 


Data Interface Options 

422A, or per User Spec 
Qualification Life (Including Particle a 
Radiation) 22.5 yrs 


Qualification Temperature 
Vibration 
Parts Program 
Reliability 


Mounting Asymmetric 3 Point Flange 
Mount 


Table 5-2: The Earth Sensor / Sun Sensor Assembly Specifications 





4. The Reaction Wheel Assembly 


a. Application 

The Type E Reaction Wheel Assembly provides a reliable source of 
reaction torque and angular momentum storage for attitude control of medium to large 
spacecraft. Accelerating or decelerating a flywheel with an integral motor provides a 
means of controlled momentum exchange with the spacecraft, which is advantageous for a 


variety of attitude control schemes. [Ref. 33] 
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Featuring high, bi-directional torque capability, a set of Type E Reaction 
Wheels is ideal for agile three-axis reaction control applications. With a momentum 
storage capacity of 50 N-m-s and an available reaction torque of up to 0.3 N-m, the Type 
E Reaction Wheel can satisfy a broad range of spacecraft requirements. [Ref. 33] 


b. Design Description 


1. Structure 

The Type E Reaction Wheel consists of a precision balanced solid flywheel 
suspended on ball bearings and driven by a brushless DC motor. The aluminum housing is 
vented to space to eliminate viscous drag torque from the flywheel due to windage. 
Removable domed covers provide access to the wheel interior for inspection with the 
flywheel installed. [Ref. 33] 


2. Suspension 

The flywheel is suspended on high capacity angular contact ball bearings. 
A low vapor pressure lubricant synthesized specifically for high vacuum spaceflight 
applications is used to lubricate the bearings, which incorporates a flight proven extreme 
pressure additive to prevent bearing wear while operating at low speeds below 50 rpm 
when an EHD film is insufficient to prevent ball to race contact. Labyrinth seals and 
barrier film coatings ensure long life operation by preventing loss of lubricant from the 
bearings due to evaporation and surface migration. A life test of the Reaction Wheel in a 
bi-directional reaction wheel application frequently crossing through zero speed has been 
operating successfully since April 1994. In addition, ITHACO is offering an On-orbit 
grease replenishment design for all their Reaction Wheels in order to reliably increase 
minimum ball bearing operating lifetime from eight to fifteen years. The newly developed 
device utilizes a Kapton heater and a simple thermostat circuit to drive a sealed grease 
reservoir to a specific temperature. The high thermal expansion coefficient of the grease 
results in the injection of a controlled amount of grease from the reservoir directly into the 
bearing. The system is unpressurized and contains no active moving parts, resulting in 
reliable, fail-safe performance. [Ref. 33, Ref. 34] 


3. Motor 
A discretely commutated, ironless armature brushless DC torque motor 
drives the flywheel in the Reaction Wheel. Maximum inertia to weight ratio is realized 


with the large diameter motor components and the unique feature of the ironless armature 


qe 


design which strategically places all of the mass of the iron and magnets on the rotor. The 
ironless armature design also eliminates cogging for smooth transitions during direction 
reversals. Hall sensors integrated directly into the armature sense the alternating flux field 
from the permanent magnets to provide commutation information for the motor driver, as 
well as a high resolution tachometer signal for speed and direction of rotation information. 
(Ref. 33] 


4. Motor Driver 

The motor driver responds proportional to an analog torque command and 
provides linear motor torque to the flywheel over the entire speed range. A digital 
interface per MIL-STD-1553 or other custom interface can be easily incorporated in the 
existing board space. Analog and digital tachometers are provided to accommodate a 
variety of user needs. Due to the high torque and current in the Reaction Wheel, the 
motor driver is packaged in a separate enclosed electronics box for flexible thermal 
management. The motor drivers are designed to minimize electrical interference, and have 
demonstrated superb electromagnetic compatibility by meeting the requirements of MIL- 
STD-461 and MIL-STD-462. [Ref. 33] 


5. Specifications 

The performance specifications for the Reaction Wheel (E - Wheel) are 
presented here in tabular and graphical form. Higher angular momentum capacities and 
higher torque versions can be made available as non-standard designs. Within the present 
Reaction Wheel case sizes, higher reaction torque can be accommodated with a rescaled 
motor driver. [Ref. 33] 
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ITHACO "E WHEEL" REACTION WHEEL ASSEMBLY 


Speed Range (rpm) 0 to + 3850 
Angular Momentum Capability 50 N-m-s (36.8 ft-lb-sec) 





120 gram-cm? (07-in?) 


Table 5-3: Ithaco Reaction Wheel Specifications 


300 mN-m (42 07-in) 
Weight of Wheel Assembly 10.5 kg (23.1 Ib) 










5. The Satellite Inertial Reference System 


a. Heritage 

Honeywell has been involved in many successful space programs including 
launch vehicles, satellites, space probes and manned space flight. Recent experiences have 
included the TOS , Centaur launch systems, MSX and Space Station Freedom. 

The Miniature Inertial Measurement Unit (MIMU) product is an 
integration of components from several Honeywell programs into a standard product. 
These other programs have provided the hardware and software knowledge and 
experience to develop a family of systems consisting of three axis rotation sensors and 
optional three axis linear acceleration sensors in a small, lightweight radiation hardened 
product. The MIMU uses the Honeywell GG 1320 ring laser gyroscope which is in high 
volume production the their Minneapolis ring laser gyro production facility. 


b. Miniature Inertial Measurement Unit 

The MIMU family was designed to be a low cost, low risk radiation 
hardened product, to be used as the inertial measurement system in a broad range of 
satellite applications. The system contains three GG 1320 ring laser gyros mounted 
orthogonally as the rotation sensors and three optional Allied Signal QA 3000 
accelerometers as the linear acceleration sensors. A block diagram of the MIMU ts shown 
in Figure 5-5. [Ref. 35] 
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Figure 5-5: The Miniature Inertial Measurement Unit 
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The Honeywell MIMU Model YG9666 B 
250 in® 
7.6 |b 
33.7 Watts 
roy S 


N 
O 


oan 
Ia 


Weight (w/o magnetic shield) 
Power 
Rate Capability (deg/sec) 
Acceleration 
Life (yrs, operation dependent) 
Gyro-GG1320 
Bias Stability 
(1o, deg/hr; short term) 
Angular Random Walk 0.005 
(deg/ Vhr ) 


Scale Factor Error 


0.008 


(lo, ppm; short term) 
Readout Noise 
(urad, Io) 


Accelerometer QA 3000 
Bias Stability 0.04 
(1o, mg; short term) 


Scale Factor Error 


(lo, ppm; short term) 
Output Data 

Angle LSB (urad) 

Velocity LSB (mm/sec) 
Temperature 

Operating (°C) -30 to +70 

Non-Operating (°C) -45 to +75 
Radiation (krads total dose) 100 

Table 5-4: The MIMU Performance Parameters 


joa 
SS 


0.0753 


The MIMU offers excellent long term stability, very high reliability and a 
wide dynamic range. Recent reliability studies of field life data has shown the life of the 
gyroscope to be in excess of 300,000 hours. This gives the MIMU a high probability of 


success for missions of 15 years or more. [Ref. 35] 


The processor electronics assembly consists of the Honeywell radiation 
hardened GVSC, RH-1750 processor, a 16 bit processor utilizing the full MIL-STD- 
1750A instruction set. The software performs the functional control as shown in the block 
diagram of Figure 5-6. This assembly also contains a standard RS422 interface operating 
at a 200 Hz data rate. MIMU can also be configured to use an optional MIL-STD-1553B 
interface. [Ref. 35] 

The standard RS422 interface consists of ten 16-bit data words. These 


data words are shown in Table V-S. 


Frame Timer (MSH) 
Frame Timer (LSH) 







X Compensated Gyro Angle 
Y Compensated Gyro 
Z Compensated Gyro Angle 






Table 5-5: MIMU Output Data Message 








The data rate is 1 Mbit/sec. The multiplexed BIT Status word consists of 
256 words multiplexed at one word output each 200 Hz period. The BIT status word 1s 
identified by the LSH of the frame timer word. [Ref. 35] 

The accelerometer electronics are also located on the processor assembly. 
The QA3000 Q-Flex accelerometer is a torque-rebalance sensor with fused quartz flexures 
with a permanent magnetic torquer, a capacitive pickoff system, and self-contained servo 
electronics. The accelerometer outputs a current proportional to the acceleration. The 
accelerometer electronics on the processor assembly digitizes this current (three channels) 
which is then read by the processor. Mathematical modeling of bias, scale factor, and axis 
misalignment is performed over the temperature range to enhance the performance. [Ref. 
35] 

The gyro electronics assembly contains the electronics for control of the 
three Ring Laser Gyros. The RLG readout post amplifier provides the gain to amplify the 


gyroscope's photodetector outputs to levels which can be digitized. These readout pulses 


are then accumulated and are read by the processor. Pulses are accumulated for positive 
angular rotation and for negative angular rotation. Net rotation is then calculated by the 
processor. [Ref. 35] 

To over come a lock-in mechanism at low rates, the gyro is dithered, which 
involves subjecting the gyroscope to a small angular oscillation through a controlled 
mechanical resonance. The MIMU gyros are oscillated at three separate frequencies 
(nominally 525, 575, and 625 Hz). The baseline means of removing the dither oscillation 
from the gyro signal is through a procedure called "dither stripping". The procedure 
involves measuring the angle between the gyro block and the gyro mounting base and 
subtracting this angle from the gyro output signal. The output of a piezoelectric sensor 
attached to the dither spring is passed through an A/D converter to the digital processor, 
scaled and corrected for phase, and subtracted from the gyro signal. The result is a wide 
bandwidth signal with the content at the dither frequency removed. An alternate means of 
removing dither is filtering. Filtering is used when low readout noise is required, and the 
corresponding latency penalty can be tolerated. The dither control electronics are also 
located on the gyro electronics assembly. [Ref. 35] 

The gyro electronics assembly also contains the path length control (PLC) 
circuitry. In order to maintain the RLG laser at its peak operating point (peak power), 
PLC is utilized by using a PZT which converts a voltage input to a path length distance. 
The lasing power is at maximum when the laser beam path length is equal to an integral 
number of wavelengths. [Ref. 35] 

The gyro electronics assembly also contains the RLG current control 
circuitry. The current in each leg of the RLG is precisely controlled to minimize 
temperature and gas flow induced bias changes. [Ref. 35] 

The Low Voltage Power Supply modules convert the spacecraft +28 Vdc 
to the isolated secondary voltages required for operations (+5 Vdc, +15 Vdc, and -15 
Vdc). [Ref. 35] 

The High Voltage Power Supply (HVPS) module provides +2200 Vdc to 
the gyroscope cathodes during startup. After startup the HVPS provides -750 Vdc to 
maintain the laser currents. The PLC circuitry requires +250 Vdc provided the HVPS. 


[Ref. 35] 


6. HD-1003 Star Tracker 
The HD-1003 Star Tracker, produced by Hughes Danbury Optical Systems, Inc. is 
designed to meet current and future needs for a lightweight, reliable and accurate star 


tracker. A simple optical system, ASIC based electronics and high throughput firmware 
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architecture make the HD-1003 mission adaptable, even while in operational orbit. [Ref. 


36] 


The HD-1003 attitude determination system provides continuous three-axis 


attitude data over the entire orbit. An autonomous all-stellar attitude determination 


system is available, minimizing spacecraft computation and ground support. Since the 


HD-1003 has no moving parts, high reliability and long life are achieved. Further, since 


the HD-1003 system can provide multiple star updates every one tenth of a second for 


LOS rates up to 1.0 deg/sec, reliance on the inertial reference unit can be reduced or 
eliminated. [Ref. 36] 
Some key features of the HD-1003 include: 


A. Performance 


lr 


. SIX arc-sec rms per star 


Six-star tracking capability 


Two arc-sec rms multi-star average 


. 10 Hz update rate 
. 8x 89 or 20 degree circular FOV 


. Parts count has been reduced by 5 times. 


B. Growth 


ie 


12 MIPS throughput with 50% margin 


2. Expandable memory 


3. 
4. 


Uplinkable code 


Small or large space-craft compatible 


C. Survivability 


le 
. Built-in radiation shielding 


A B&W N 


All electronic parts are hardened to 100 Krad or higher 


. Space qualified proton/debris rejection software 


. Expected life on orbit is 15 to 18 years 


HDOS has completed a 2 year test program to verity the HD-1003's 


performance in high radiation environments. 


6. Test results and design have been independently verified by Lockheed 


Martin, Aerospace and the Naval Research Lab. 


eo 


D. Producibility 
1. Modular design 
2. ASIC-based electronics 
3. No moving parts, sun safe 


4. Safety engineered 


The HD-1003 Star Tracker benefits from over Twenty years of experience in 
leading attitude determination technology. Hughes Danbury produced the world's first 
CCD-based star sensor to be qualified in space, the first "multi-star" tracker, and the first 
CCD upgrade to NASA's Standard Fixed Head Star Tracker. The Hughes Danbury team 
continues to satisfy evolving demands by utilizing the latest technology and innovative 
packaging. [Ref. 36] 


HOUSEKEEPING | [TIMING CONTROLLER 
DATA ACQUISITION] | AND CCD CLOCKING ee 


SYSTEM CPU 
PROGRAM AND DATA 
MEMORY 


[\ ee SIGNAL 
E 2 PROCESSOR 





PIXEL DATA 
AND ADDRESS 
PROCESSOR 














\ 


SECONDARY POWER 
TO ALL LOADS 


COMMAND / DATA POWER 
INTERFACE SUPPLY 
HOST INTERFACE 


Figure 5-6: CCD Focal Plane and Microprocessor Controller Sample and Centroid up to 


Six Star Images Simultaneously. 
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Size (Win iene ehed en in een 7e6O x 44 in 
Weight (with light shade) 7.1 Ib 

Power (average at 28 Vdc) 9 Watts 
Communication Interface MIL-STD-1553B 
Field of View (deg) 20° circular 


Performance Category HD-1003 Narrow HD-1003 Wide FOV 
FOV 


| Sensitivity (5000 K star) (m, ) +6.0 | 44.7 
} Overall Accuracy, Each Axis 

* Pitch / Yaw (per star) (arc-sec, 6 30 
| rms) 2 10 

* Multi-star average (arc-sec, rms) 40 50 

* Roll (5-stars) (arc-sec, rms) 25 0.25 


* Magnitude (+) 
Update (frame) Rate (Hz) 


Stars Simultaneously Tracked 6. ee 
3 
2 
2 


4 i 
cs 
Bright Object Rejection Angle** 
* Sun (deg) a 
* Earth (deg) Se 
* Moon (deg) Se 


4 

40° 

30° 

255 
Acquisition Time (6 stars) (sec. 1-sigma) i; ii 
Mean Time Between Failures (Hours) 


Environments: 





Temperature 


* Nominal ( °C) -5 to +45 -5 to +45 
* Survival ( °C) -15 to +55 -15 to +55 
Vibration (Random) (g rms) Wee We 
EMC Level MIL-STD-461C MIL-STD-461C 
** Light shade is modular to accommodate various mission requirements 
Table 5-6: HD-1003 Star Tracker Performance Capability 





7. Global Positioning System 


A Global Positioning System receiver is required to provide the spacecraft with an 
autonomous capability to accurately determine its position. The Motorola GPS 
Spaceborne Receiver system was selected for this mission. This system is capable of 
determining the spacecraft's position to within 15 meters for Autonomous Position and 2 
to 5 meters for Differential Position; provided that the Selective Availability is not 
implemented. This system also provides the ability to determine spacecraft attitude to 
accuracies of 0.1° in pitch, roll and yaw. In addition, the system can determine pitch, roll 
and yaw rates to accuracies of 0.5°/second. These accuracies are not sufficient to meet 
the crosslink pointing knowledge requirements but are suitable to provide a backup source 
of attitude information for basic spacecraft control functions in the event of failures in the 
primary. Table 5-7 identifies the characteristics of the Motorola GPS Spaceborne 
Receiver. 

Additional key features of the Motorola GPS Spaceborne Receiver include: 


Space qualified digital design. 
. Full spaceborne doppler capability. 
. Autonomous operation with doppler capability at 1 second rate. 


. Non-destruct doppler measurements at 1 second rate. 


nA fF WN = 


. Radiation hardened 1 PPS clock output synchronized to GPS time. 
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Multiple Configurations Available 
* L1, C/A and P-code 

* 1.1/L2, C/A and P-code 

* Y-code-Anti-Spoof and Selective Availability 
* Multiple Antennas 

* 6 to 18 channels 








Physical / Environment 


Size Gl x 5.2) ee 
Weight 3.7 Ib 
DC Power 


| Vibration 25 grms 
Shock 3000 g at 10 kHz 
| Temperature 20m CiCan ae 


Performance Characteristics 
Receiver Accuracies (over C/N | Range of 32 to 42 dB/Hz): 
Pseudo-Range (10 sec) 0.2 m 
Pseudo-Range Interchannel Bias <0.1m 
Carrier Phase (1 0.01 m 
Navigation Solution Accuracy (GDOP = 3, Selective Availability not 





























sec) 


implemented): 
Autonomous Position 
Differential Position 
Time Offset 

Time-to-First-Fix 


Mean TTFF 
Maximum TTFF (P fo 0.9) 


<15m(lo) 
2 to 5 m(lo) 
100 ns 


< 60 sec 

< 90 sec 

< 10 km, 10 m/s 
<1x10” 


User S/C Ephemeris 
User Clock Error 
Orbital Dynamics (typical / spacecraft orbit): 


Velocit 8,000 m/s 
Acceleration 


Table 5-7: Characteristics of the Motorola GPS Spaceborne Receiver 
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8. Sun Sensor 


a. Overview 

Two Sun sensors will be used; one on each solar array substrate assembly. 
These two sun sensors will be the primary method for determining the yaw angle and solar 
array rotation angle required to maintain solar tracking. The Sun Sensor Assembly 
selected is a two axis digital sun sensor which offers unprecedented performance in terms 
of accuracy, mass, power consumption, reliability and automation. The sensor 
autonomously determines two-axis attitude information with respect to the sun line from 


low Earth through super geosynchronous orbits. 


b. Theory of Operation 

The Sun Sensor Assembly incorporates two static, mutually orthogonal 
sensing (axes) arrays to detect the orientation of the sun. Each axis is comprised of a slit 
allowing a plane of light to pass through it intersecting with an orthogonally mounted, 
high resolution, linear photodiode array. The point of intersection, which indicates sun 
angle, is obtained by reading the energy stored in the array. 

The Sun Sensor Assembly typically operates with sun locations of + 64° in 
azimuth and + 64° in elevation. All processing and calibration data utilized to derive sun 
azimuth and elevation from the Sun Sensor Assembly occurs on board the sensor. Output 
occurs over a digital interface governed by the MIL-STD-1553 data bus, EIA RS-422A, 
or user specified interface. Upon asserting a data request the SSA outputs all SSA 
information including diagnostics and status information. Units have been provided for 
the MSTI, APEX, and Sea Star Programs. 
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VI. YAW CONTROLLING 


On orbit pointing control utilizes three reaction wheels to provide roll, pitch and 
yaw attitude control. An earth sensor provides two axis attitude data for roll and pitch 
and a sun sensor provides one axis attitude data for yaw. The yaw control is similar to 
both roll and pitch except that it must utilize additional information in order to maintain 
solar array sun tracking. Also, changes is yaw do not adversely effect the communications 
payload as long as pointing error limits about the roll and pitch axis are maintained. In this 


section the peculiarities of yaw steering will be discussed. 


A. YAW DISTURBANCES 

Due to secular disturbances, momentum accumulates in the wheel, and wheel 
speeds will build up with time. In normal operations these secular torques in the yaw 
direction will be compensated by operation of the magnetic torque rods In the back-up 


mode momentum would be dumped through a series of minimum impulse thruster pulses. 


B. SUN AND NADIR POINTING GEOMETRY 

Figure 6-1 illustrates the Sun and Nadir pointing geometry for the MUS 
spacecraft. represents the azimuth and 6 the elevation of the sun line relative to the 
spacecraft orbit plane. B is referenced to the X axis. The X - axis coincides with the 
velocity vector in the local geodetic reference frame (i.e. orbit nominal frame). The 
commanded yaw and array angles vary as a function of f and 6 due to the effects of the 
space vehicles orbital motion, precession of the orbit plane, and movement of the sun in 
the earth centered inertial frame. These commands maintain the normal to the solar panels 
directed toward the sun and the sun in the spacecraft roll/yaw plane. Except for very low 
sun angles, when 6 < 5.0°, the -X side of the spacecraft is always oriented towards the 
sun. [Ref. 27: p. 458] 
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f, 


Figure 6-1: Pointing Geometry 


C. SUN AND NADIR YAW POINTING STRATEGY 

The Attitude Dynamics and Control! Subsystem is designed to track the nominal 
sun/nadir yaw trajectory for sun angles (6) greater than 5.0° from the orbit plane. The 
closed loop controller uses sun sensor data as the yaw attitude reference. When the 
spacecraft is in the sun/earth co-linearity region, defined when 6 is less than 15° and the 
projection of the sun line in the orbit plane is within 15° of the nadir vector, yaw sensing 1s 
restricted due to either poor geometry or an eclipse condition. In these regions, a yaw 
momentum tracking tach loop is used to track an ideal yaw momentum profile. Predicted 
Euler coupling torques are fed forward to achieve a high degree of accuracy, minimizing 
yaw pointing error at exit from the region. [Ref. 27: p. 460] 
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When the sun is within 5.0°, nominal, of the orbit plane, the ADCS enters a 
constant yaw mode which maintains the pitch (array rotational) axis on the orbit normal 
and an orbit rate rotation about the pitch axis. The constant yaw mode minimizes 
Reaction Wheel torque and momentum demand while maintaining specified yaw pointing. 


Table 6-1 summarizes the sun and nadir yaw pointing strategy. [Ref. 27: p. 460] 


Sun Rotation Sun to Yaw Control Remarks 
Tracking | Angle B(deg) | Orbit Plane Reference 
Angle (6) 
All B angles [Olea ST/Yaw SS Closed loop yaw control using 
Output Null Star Trackers (ST) / SS derived 





Position yaw & yaw rate 


Ideal IB - 90°| = 15° 15° > |d| = ST/Yaw SS Closed loop yaw control using 
Yaw and 5.0° output null ST/SS derived yaw & yaw rate 
Tracking | |B - 270°| = 15° position 
Mode 
IB - 90°] << 15° | 15°>|d}2 | ST/Calculated | Poor yaw signal from SS when 
or 5.0° yaw cmd, no | {8} < 15° due to array position or 
|B - 270°| < 15° SS yaw no SS output due to eclipse 
sensing 
Ip - 90°] = 15° [5] < 5.0° ST/Yaw SS _ | Constant pitch rate control; 
and output biased | Closed loop yaw control using 
|B - 270°| = 15° at 6 ST/SS derived yaw and yaw rate. 
Constant 
yaw |B - 90°| < 15° |5| < 5.0° ST/Zero yaw | Constant pitch rate control; Feed 


tracking or angle; no SS__| forward tach loop yaw rate 


mode |B- 270°| < 15° yaw sensing | tracking control; Poor SS yaw 
signal when |B| < 15° due to 
array position or no SS output 
due to eclipse; Minimum yaw 


control torque required. 





Table 6-1: Sun and Nadir Yaw Pointing Strategy 
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D. SUN AND NADIR ARRAY CONTROL STRATEGY 

The ADCS is designed to track the nominal sun/nadir array trajectory for sun 
angles (6) greater than 5.0° from the orbit plane. The closed loop controller uses course 
sun sensor data as the array loop null reference. When the spacecraft is in eclipse, a 
calculated nominal array trajectory drives the array loop with position feedback provided 
by the measured array angle. When the sun is within 5.0°, nominal, of the orbit plane, the 
ADCS enters a constant yaw mode. The array control strategy is unchanged, but the 
constant yaw geometry results in an orbital rate rotation of the array. The transition into 
and out of the constant yaw mode is limited to the segment of the orbit with sun rotation 
angles between 120° and 240°; 1e. between 2pm and 10pm spacecraft local time, 
precluding the need for yaw inversion maneuvers. Table 6-2 describes the sun/nadir array 
control strategy. [Ref. 27: p. 462] 


Solar Array | Eclipse Sun to Orbit Array Control 
Tracking | Condition | Plane angle 6 Reference 
Mode 
No |S] = 5.0° ST/Pitch SS output | No eclipse will occur 
Ideal Sun- eclipse null position when |d| > 15° 
Nadir 
Array The error between 
Tracking Eclipse lise 22 (6) ea 5008 Calculated ideal measured array 
array angle command position and 
command reference 
is used for control 
No |5| < 5.0° ST/Pitch SS output Array sun pointing 
Constant eclipse null position error 6 
Array 
Rate |5| < 5.0° Calculated array Measured array 


Tracking Eclipse angle command for position used for 


constant yaw control; sun pointing 





tracking error equals 6 
Table 6-2: Sun and Nadir Array Control Strategy 
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E. NORMAL, ON-STATION CONTROL FUNCTIONAL BLOCK 
DIAGRAM 

The Attitude Dynamics and Control Subsystem provides fully autonomous three 
axis attitude determination and control throughout the Space Vehicle's operational 
lifetime. Three of four redundant reaction wheels control body attitude. and redundant 
magnetic torquing coils control system momentum. There is no thruster utilization for the 
Normal, On-Station mode of operation. Thrusters are used primarily to provide A V and 
control actuation during stationkeeping maneuvers, but also provide backup momentum 
unloading in the event of a failure of the magnetic system. Also, control loop 
compensation is provided in the attitude control software for Euler coupling torques, 
magnetic and thruster desaturation torques, and wheel friction torques to optimize 
pointing and reduce jitter. [Ref. 27: p. 266] 

While the system is fully autonomous, the Control Segment has full override 
access to the ADCS equipment configuration and operating mode. The internal state of 
all software flags is available in telemetry and any control mode or functions may be either 
inhibited or forced by a ground override of the appropriate flag(s). Any control algorithm 
may also be reprogrammed by ground command. The system is capable of operation for 
periods of at least 210 days with no ground intervention. Figure 6-2 illustrates the 
Normal, On-Station Functional Block Diagram. [Ref. 27: p. 466] 
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Figure 6-2: Normal, On-Station Functional Block Diagram 


F. ATTITUDE DETERMINATION SOFTWARE FUNCTIONAL 
OVERVIEW 


The attitude determination software comprises two primary functions; namely, 
attitude sensor data processing and payload ephemeris data processing. The attitude 
sensor processing conditions the raw data input from the various attitude sensors, 
processes the data to derive three-axis spacecraft attitude and attitude rate, and generates 
the yaw and solar array drive reference commands to drive the respective control loops in 
the co-linearity regions and for low sun angles (6<5.0°). The ephemeris processing, 
which uses Navigation Payload data as the primary reference and ground uplink data as 
the backup reference, generates sun angle data for sun/nadir pointing mode operation, 
sun/moon interference predictions for earth sensor processing, and a magnetic field model 
for the momentum unloading software. Figure 6-3 is a Functional Overview of the 
Attitude Determination Software. [Ref. 27: p. 468] 
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G. SUN ANGLE ESTIMATION (8B AND 5) 
The ADCS processes both Navigation payload data as well as ground uploaded 


ephemeris to provide a reliable basis for sun/nadir pointing. The spacecraft ephemeris 
obtained from the on-board navigation payload and the sun ephemeris generated on-board 
are used to compute the sun rotation (azimuth) angle (B) and the out-of orbit plane sun 
angle (6) from which eclipses times and yaw/array control phase is determined. [Ref. 27: 
p. 470] 

The backup mode determines 6 using the orbit right ascension of ascending node, 
Q. and the orbit inclination, 1, along with the on-board generated sun ephemeris. B is 
estimated from the calculated 5 angle and the array position when the array position 
relative to the yaw axis is suitable and attitude pointing is within tolerance. Yearly 
variations of (2 and 1 are represented by an on-board model, thus each upload covers a 
time span of one year. To satisfy the autonomy requirements, however, uploads are 
required every few months. The backup estimator is run at all times, but the estimated B 
and 6 are only used when payload data is unavailable or invalid. Figure 6-4 is the flow 
diagram for the estimation of Sun angles B and 6. [Ref. 27: p. 470] 
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Figure 6-4: Sun Angle Estimation 
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H. YAW DATA PROCESSING 


The ADCS utilizes a two-axis sun sensor for yaw attitude control and solar array 
drive control. The Sun sensor is comprised of two sensing units, one located on each 
solar array, which provide redundant detector sets and processing electronics for each 
sensing axis. Data from the two detector units is differenced to provide a measurement 
that is linear on each axis over approximately + 40° and covers a full range of + 125°, 
minimum, per axis. For yaw determination, the sensor data must be weighted as a 
function of the solar array orientation. [Ref. 27: p. 472] 


10] 


This Page Intentionally Left Blank 


Vil. SIMULATION OF THE ATTITUDE DYNAMICS AND 
CONTROL SUBSYSTEM 


A. PROBLEM STATEMENT 
The MUS spacecraft is a Medium Earth Orbit satellite which has the following 


moments of inertia of: 


me Sree nt 
Ue = 14,083.6 kg . mr 
Dyaw = 14,113.7 ke «me? 


The pointing accuracy for the spacecraft is: 


a) 0.1° for Roll and Pitch 
b) 1.0° for Yaw 


The disturbance torques on the spacecraft are primarily due to solar pressure, and 


can be approximated by: 


T. = T.,, = (1x 10”) *cos(@,t) N-m 

ie ns 10m cos(yt) Ni 
= 

T, = T,,, = (1x 10~) * cos(@,t) N-m 


The thruster misalignment torque about any axis can be approximated by: 


en N-m 


Finally, the design of the attitude control system for the spacecraft must meet the 


given specifications when the environmental torques are imposed on the spacecraft. 
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B. THREE AXIS STABILIZATION 


The control torques along the axes of three-axis-stabilization systems are provided 
by various combinations of momentum wheels, reaction wheels, and thrusters. Broadly, 
however, there are two types of three-axis-stabilization systems: a momentum biased 
system with a momentum wheel along the pitch axis, and a zero-momentum system with a 
reaction wheel along each axis. [Ref. 25: p. 130] 

In momentum wheel systems, the angular momentum along the pitch axis provides 
gyroscopic stiffness. In these systems, the pitch and roll axes are controlled directly and 
the yaw axis is controlled indirectly due to a yaw sensor. In reaction wheel systems, all 
three axes are controlled independently, thus requiring a yaw sensor also. The momentum 
wheel can be fixed or gimbaled about one or two axes. The control torque along the pitch 
axis 1s provided by the change in the speed of the momentum wheel. The torque along the 
roll axis is provided by thrusters, a reaction wheel, or by changing the gimbal angles in the 


gimbaled momentum wheel. [Ref. 25: p. 130] 


1. Small Angle Approximation 

When the satellite's body frame orientation remains very close to the orbit nominal 
frame, a small angle approximation is achievable as well as being able to disregard higher 
order terms. This greatly simplifies the problem. A linearized analysis of the three-axis- 
stabilization system is given in this section. The linearized equations of motion of a three- 
axis-stabilized spacecraft are derived first. The angular velocities of the spacecraft can be 


expressed in terms of the orbital rate, @, , and the attitude error angles, yw, 0, and @, 


which are known as yaw, pitch, and roll errors, respectively, as follows: 
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x | 0 -sin(@) |/¢ 
Dia soyi—1\0 cos(g) cos(A)sin(¢) |; + — 
@ 0 -sin(g) — cos(A)cos(¢g) || v 


cos(9) sin(y) d-aow 
@,\co(g)cos(y) + sin(¢)sin(6)sin(y) +} =40-a, (7-1) 
-sin(¢)cos(y) + cos(¢)sin(4)sin(y) yt+o$ 


The angular momentum of the system can be written in the form 


ee eH, (7-2) 


w 


where, 


H,=[h, h, he] 
e r (7-3) 
H,=(1,0, 1,0, 1.0, | 


J 


are the angular momentum of the wheels and of the spacecraft, respectively, where 
I.,1,,,and I,are the principal moments of inertia of the spacecraft and 
@,,@,,and w, are the angular velocity components of axes x, y, and z as given by Eq. 


(7-1). Introducing Eq. (7-2), in conjunction with Eqs. (7-3) and (7-1) into 


—~ 


M. = Herei+ @ x H, and ignoring nonlinear terms, we obtain the moment equation in the 


form 


h, + 1,.(¢-@,w)+h(0-0,)-hv + 0,9)-(1. - 1, )o.(v + 0.9) 
Mz h, - 1,0 a h,(w + @,¢)- h(¢ = ow) (7-4) 
h, + LAW aE 09) t h,(¢ = ow) = h(@- @,) = (7,, ele jo,(¢ 5 ow) 
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2. Large Angle Requirements 

Unfortunately, the yaw steering requirements of the satellite does not allow us to 
either make small angle approximations or to disregard higher order terms. This greatly 
complicates the equations of motion which are presented below. As previously discussed, 
the angular velocities of the spacecraft can be expressed in terms of the orbital rate, @, , 


and the attitude angles, ¢, 8, andy, which are known as roll, pitch, and yaw respectively, 


as follows: 
Oo, I 0 — sin(0) p 
O=4@,7=|0 cos( ¢) cos(@)sin(¢) Q\— 
@ 0 -sin(¢) — cos(A)cos(¢) |] ¥ 


cos(0)sin(y) 
@,, co ¢)cos(w) + sin(¢)sin(6)sin(y) 
—sin(¢)cos(y) + cos(¢)sin(@)sin(y) 
¢— wsin(@) — w, cos(9) sin(y) 


= 4 Acos(¢) + wcos(@)sin(¢) — w, cos(¢)cos(w) — a, sin(¢) sin(9) sin(y) (7-5) 
—Osin(¢) + wcos(@)cos(¢) + w, sin(¢)cos(y) — @, cos(¢)sin(@) sin(y) 


The angular momentum of the system can again be written as stated in Eq. (7-2) 
H=H, +H, (7-2) 


where, 


cS} 
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are the angular momentum of the wheels and of the spacecraft, respectively, where 


I..1,,,and I, are the principal moments of inertia of the spacecraft and 
@,,@,,and w@, are the angular velocity components of axes x, y, and z as given by Eq. 


(7-5). Introducing Eq. (7-2), in conjunction with Eqs. (7-3) and (7-5) into 


M, = Herei+ @ x H. and keeping all terms, we obtain the moment equations in the form: 


h, + 16+ Ocos(dyh, + dcos(¢) lw + ywcos(@)sin( g)h, + w* cos(9) sin(¢) I. — 
@, cos(g)cos(y)h, — @, cos(¢) cos(y)/_y — a, sin(¢) sin(6) sin(y )h, — 
@, sin(g)sin(@) sin(y)/_w + Osin(g)h, + & sin(¢)1 w ~Wcos(@)cos(@)h, — 
y cos(9)cos(¢)/,,0- a, sin(¢)cos(y)h, — w, sin(¢)cos(y)1,,.0 
+ w, cos(g)sin(8) sin(y)h, + @, cos(g)sin() sin(y)/,,0 


h, + 1,,0—Osin(g)h, — Osin(¢)1,.¢ + cos(8) cos(¢)h, + y cos(P) cos(¢)1,.4 
= + @,sin(¢)cos(y)h, + w, sin(¢)cos(y)/,.¢ — a, cos(¢)sin() sin(y)h, — 
wo, cos(¢) sin(9)sin(y)/,,.¢— gh, - ¢l..w + wsin(O)h, + yw’ sin(O)I_ + 
@,cos(@)sin(yw)h, + @, cos(@)sin(y) ly 


(7—6) 


h, + lw + gh, + $l,,0- ysin(@)h, — ysin(@)1,,0— @, cos(A)sin(w)h, - 
o,cos(@)sin(y)1,,0 — Ocos(¢)h, — Ocos(¢)1,,.¢ — wcos(6)sin(g)h, — 
yw cos()sin( ¢)/,.¢ + @, cos(¢) cos(y)h, + w, cos(¢)cos(y) 1, + 
@, sin(¢)sin(O)sin(y)h, + @, sin(¢) sin(9) sin(y)1,,¢ 


The external moments arise from three major sources: the gravitational gradient, 
the solar radiation pressure, and the control moments from actuators. Denoting the 


individual torques by M., M.. and M., respectively, we can write 


M=M,+M,+M.. Gey) 
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where. 


d(1..- 1} 
Mz, = 3035 AL. — 12) (7-8) 
0 


Different types of active control systems have different key parameters and 
algorithms. Frequently, three-axis control can be decoupled into three independent axes. 
The most basic design parameter in each axis is its position gain, K,. This is the amount 
of control torque which results from a unit attitude error and can be expressed in Nm/deg 
or Nnvrad. The position gain is selected by the designer and must be high enough to 
provide the required attitude control accuracy in the presence of disturbances, or 


= li 0.” where K,, is position gain, 7/, is peak disturbance torque, and @, is allowable 


attitude error. In its simplest form, a spacecraft attitude control system can be represented 


in the s-domain as a vA 2 plant and may be controlled by a proportional plus derivative 


(PD) controller where 7. = K,0, + K,6.. The position gain, K,,controls system 
bandwidth and the rate gain, K,, controls damping. [Ref. 6: p. 364] 


The value of the position gain also determines the attitude control system 
bandwidth and speed-of-response. The bandwidth is given by , = (K, ia *, where / is 


the spacecraft moment of inertia. The bandwidth defines the frequency at which control 
authority begins to diminish. Attitude control and disturbance rejection are effective from 
0 frequency (d.c.) up to the bandwidth. Speed of response is approximately the reciprocal 
of bandwidth. Note that position gain is inversely proportional to allowable error and 
bandwidth is proportional to the square root of position gain. Therefore, high accuracy 
implies high position gain and high bandwidth. However, high bandwidth may cause 


bending resonances to affect control system performance. [Ref. 6: p. 365] 
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With the relations given, the system designer can estimate required position gain 
from his estimates of disturbance torque and accuracy requirements. He can use this 
estimate to compute control system bandwidth. This allows him to specify minimum 
bending frequencies as discussed below. [Ref. 6: p. 365] 

In defining algorithms for the control system, we must also consider whether the 
vehicle will have flexible-body effects that can make the vehicle unstable. Spacecraft with 
flexible appendages such as antennas, booms, and solar panels may produce slight warping 
at their natural frequencies. Control torques and external-disturbance torques will cause 
structural vibrations, in some cases close to or within the control system's bandwidth. The 
lowest natural frequencies of flexible components should be at least an order of magnitude 


greater than the rigid body frequencies before we can neglect flexibility. 


C. THREE-AXIS REACTION WHEEL SYSTEM 


A three-axis reaction wheel system can be considered as a combination of three 
independent roll, pitch, and yaw control systems, see Figure 7-1. Each axis is controlled 
by varying the speed of the reaction wheel in response to the attitude error. The system 
requires a reaction wheel with zero nominal angular momentum and an attitude sensor for 


each axis. [Ref. 25: p. 149] 
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Figure 7-1: Block diagram of a three-axis reaction control system 


1. Small Angle Approximation 


In this section the control of the spacecraft is explored when small angle 
approximations are made which would be the case if yaw steering was not being 


performed. 
Introducing Eqs. (7-4) and (7-8) into Eq. (7-7) and neglecting the transverse 


inertia of the wheel, the equations of motion for a three-axis reaction wheel system are 


M+ M,, =1,,¢+ [4e;( 1, -I_)- o,h,|9 + \-h, ~ (Tee ~ Ly + 1.) |W 
+ Oh. —@,h, + h, 


M,,+ M, =1,,0+ 30301. — 12) + @,h.b — ht ohy + hw, +h, (7-9) 


M+ Me = 1a +105 (Ip ~ Lee) - Osh, |v +[h, + @,(Tee + Le Iy)|d 
—~h.O+a,h, +h. 
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Although Eqs. (7-9) are coupled, because h,,h,,h,,and @, are small, the 


coupling terms are small. If the coupling terms are neglected, the equations of motion 


about the pitch, roll, and yaw axes become independent, and hence they can be controlled 


independently. Therefore, the equations of motion become: 


M,,+ M,, = 1,.6+h, = T, 
M,+M, =1,0+4, =T, 
M.+M, = Ly +h, z ib 


(7-10) 


Then, the control torques are applied by letting the rates of change of the angular 
moment of the reaction wheels have 


h, = K,(7,0 + é) 
h, = K,(r,0+ 6) 


(7-11) 
1, = K,(t,w+y) 
Making the appropriate substitution leads to the following: 
T, = 1,.0+K,| 76+ 9] 
T, = 1,0+ K,|7,0+ 6] (7-12) 


T, = 1, +K[t.y+y] 


Next, entering the Laplace domain yields: 


11] 


T.(s) 


As) = [78° +K,7,5+K, | 
Hs) = ee (7-13) 
[yy5° + Kt + K,| 


r(9) 


v4s) = eg + Kz,s+ K,| 


The next step is to determine the position gains and the time constants. For the 


roll axis we can say: 





¢.= lo, (7-14) 


Then, using the final value theorem, ¢,. = Lims ((s)] , we can state 


(2x 10°) 
03. = 
K, 
and since, ¢,. < 0.001745 rad = K,>0.1146. Therefore, we can let 
Kae 





Next, since w, = - = 0.0061227 rad/sec=9.7446e-4 hz, 


XX 





and, since 7, = 261 =, tT, = 2060. 
Oo 


r 


In a similar manner the remaining values were determined and are listed in Table 7-1. 


Pe 











ROLL PITCH YA 
a, = 0.0061 rad/sec @, = 0.0046 rad/sec @,, = 0.00266 rad/sec 
é 


Table 7-1: Control Law Values 


lis 








| 


= | 
,= 751 


J 





The final step is to code in Matlab and simulink. The Laplace equations I used for 


this simulation are as follows: 


T,(s) - K|z,s + 1]¢(s) + | Te ae 1.,|s y(s) = Lec5” + 4a2(,, —If, }}4(s) 


T,(s)- Kea + 1]6(s) = fies +307 (1, - 1, )o(s) (7 —15) 
T,(s) _ K,|z,s + 1y(s) - 0,( Te le 1,,)s@(s) = [es + o2(1,, — 1.) }W(s) 


For the computer simulation I used Matlab and Simulink. I considered three 
situations. First, a step input due to thruster misalignment. Second, the sinusoidal input 
due to solar pressure torques; and third, an impulse input of approximately one order of 
magnitude larger than the others; see simulink diagrams and Matlab code in Appendix F. 

The first simulation was the pitch response due to a step input. As can be seen 
from Figure 7-2, a step input occurred at time, t=100. Then, the control system was able 
to counterbalance the input with in 1300 seconds and remain with in a pointing tolerance 
of 0.1°. But, this requires a continual supply of angular momentum from the reaction 
wheel. A momentum dump will need to occur after 850,000 seconds or every 9.8 days for 
the momentum wheels selected. Next, with a sinusoidal torque input the system remained 
within pointing accuracy and fluctuated between £0.05°. Finally, with an impulse input of 
approximately one order of magnitude larger than the other torques the system stabled out 
after 1600 seconds and also remained within the required pointing accuracy, see Figure 7- 


4. 
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The second simulation was of the roll and yaw components, which are coupled 
(see Appendix F). Other than the coupling difference the simulation process was the same 
as for the pitch. Graphs of the simulation can be seen in Figures 7-5 through 7-8. For the 
sinusoidal input for roll the angular momentum did not increase with time but instead was 
sinusoidal. This is significant because a momentum dump will not be required. For the 
sinusoidal input for yaw a momentum dump must occur every 800,000 seconds or every 
oe) days. 

The final simulation was that for the contingency mode of operation which 
simulated a yaw control failure (see Appendix F). For this I just used a fourth reaction 
wheel canted between the roll and yaw directions. Then, I ran all three inputs to show 
that the system was able to control the spacecraft and remain within pointing accuracies. 
This information is displayed in Figures 7-9 through 7-12. 

Finally, in estimating the fuel required for momentum dumping I tended to be 
conservative. For the estimation I am assuming that both pitch and yaw angular 
momentum must be dumped every 10.5 days and 9.25 days respectively. Therefore, we 
will need approximately 80 pulses in a year. 

Next, from the following relationships we can determine the amount of propellant 


mass used in each burn: 


Poa 
q16 
Ft 
M,=—— 
bese 


Therefore, the total propellant mass used in one year for thruster only momentum 


dumping is estimated to be 1.5 kg. 
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Figure 7-3: Pitch Response Due to a Sinusoidal Disturbance Torque. 
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Figure 7-4: Pitch Response Due to an Impulse Disturbance Torque. 
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Figure 7-5: Roll Response Due to a Sinusoidal Disturbance Torque. 
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Yaw Angular Response 


6 Yaw Torque vs Time 


Torque (N-m) 











Time (sec) 
Yaw angle psi mae 
0 1 2 3 4 5 6 7 8 9 10 

& Time (sec) 10 
3 Angular Momentum vs Time 
o 
E 
3 
= 
i 
‘x 
5 
D> 
c 
<x 

Time (sec) x 10° 


Figure 7-6: Yaw Response Due to a Sinusoidal Disturbance Torque. 
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Figure 7-7: Roll Response Due to an Impulse Disturbance Torque. 
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Figure 7-8: Yaw Response Due to an Impulse Disturbance Torque. 
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Figure 7-9: Roll Response Due to a Step Input Disturbance Torque - Back Up System. 
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Figure 7-10: Yaw Response Due to a Step Input Disturbance Torque - Back Up System. 
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Figure 7-11: Roll Response Due to an Impulse Disturbance Torque - Back Up System. 
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Figure 7-12: Yaw Response Due to an Impulse Disturbance Torque - Back Up System. 


120 


2. Yaw Steering Requirements 


In this section the attitude dynamics and control of the spacecraft is explored when 
small angle approximations are not made and all higher order terms are included in the 
equations of motion. This is a complexity that is a result of the requirements for yaw 
steering. 

Introducing Eqs. (7-6) and (7-8) into Eq. (7-7), the equations of motion for a three 


axis reaction wheel system are 


M., + M,, =1,.¢+ Ocos(¢)I_.w +  cos(6) sin(¢)I_ — a, cos(¢) cos(w)I_w — 
a, sin(¢) sin(@) sin(w) ly + e’ sin(¢) I, — wcos(@) cos(¢) 10 — 
@, sin(g)cos(y) 1,0 + @, cos(¢) sin(@) sin(y) ie - 
37, sin(¢) cos(g)[cos(O)] I. + 3«2 sin(g) cos(g)[cos(A)] I, + A, 
+ Acos(¢)h, + ycos(6) sin(g)h, — w, cos(g) cos(w)h, — 
a, sin(g) sin(@) sin(w)h, + Asin(g)h, —ycos(@) cos(¢)h, — 
w@, sin(g)cos(y)h, + @, cos(¢) sin(@) sin(y )h, 


M,, + M,, = ae ~ Osin(¢)1,,.¢ + wcos(@) cos(¢)1,,.¢ + a, sin(¢) cos(y)1,,¢ — 

— @, cos(¢g) sin() sin(y)/,.¢-d@L.wt+w sin(@)/I_. + 

@,,cos(@) sin(w)I_w — 3a@° sin(9) cos(@) cos(¢) I, + 

3@° sin(@) cos(@) cos(¢) 1. + h, ~ Osin(g)h, + cos()cos(¢)h, + 

a , sin(¢)cos(y)h, — @, cos(¢) sin(A) sin(w)h, — dh, + wsin(O)h, + 

a, cos(@) sin(y )h, 

(7-17) 

M, + M, =I, + 41,0 - wsin(6)1,,0 — @, cos(@)sin(y)1,,0 — Ocos()1,.9 — 

y cos(@) sin(¢)/,,¢ + @, cos(¢)cos(w)/,,@ + a, sin(g) sin(8) sin(y) 1.9 

— 3@? sin(9) cos(@) sin(¢)/,, + 3a? sin(A) cos() sin(g)/,, + h, + oh, — 

y sin(@)h, — w, cos(A) sin(y)h, — Ocos(¢)h, — ycos() sin(¢)h, + 

@,,cos(¢) cos(w)h, + @, sin(¢) sin(8) sin(y )h, 

These equations of motion are somewhat overwhelming. In the process of 


developing a control law the first step was to transform these equations of motion into a 


state space formulation and develop a Matlab/Simulink model for the motion of the 


spacecraft with no control. The state space formulation used was the following: 


x,=¢ 
pate 
x,=0 
x,=%,=0 
X35 = 
X,=X,=W 
X, =X, 


_=h=o= F*ecos LW + 7 cos(O)sin(g)1_, — w, cos(g) cos(y)1..¥/ — 


X,=X, 


XX 


,, sin(¢)sin(@) sin(w)/_w + @ sin(¢) I wy ~ ycos(A)cos(¢)1 ww ~ 
a, sin(¢)cos(w)1 vO + w,,cos(¢)sin(@) sin(y) 1 Ww — 

30° sin(¢) cos(¢)[cos()]° 1 +3° sin(¢) cos(¢)[cos(6)] i 

+h, + @cos(¢)h, + weos(4)sin(¢)h, — w, cos(¢) cos(w)h, — 

o, sin(¢)sin(@) sin(y)h, + Asin(d)h, — y cos(@)cos(¢)h, — 

w, sin(¢)cos(y)h, +@,cos(g)sin(@) sin(y)h, — M,,) 


(7-18) 


X, =X; = = Ff dain 1d + wcos(@) cos(¢)1,.¢ + @, sin(¢)cos(y) 1d = 


w,,cos(¢) sin(@) sin(y)/,.¢—-dl_.W + @, cos(@)sin(w)1_w — dh, 

— 3«;, sin() cos(9) cos() I_, + 3a, sin(A) cos(A) cos(¢)/,, +h, 

— Osin(¢)h, + wcos(@)cos(¢)h, + a, sin(d)cos(y)h, + wy’ sin(A)I_ — 
@,,cos(¢)sin(@) sin(w)h, + ysin(@)h, + @, cos(@) sin(y)h, — M,,,) 
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=i. = = @i,0 ~ wy sin(9)/,,0 — @, cos(9)sin(y)/,, 8 — Ocos(¢)I,.¢ 
y cos(9)sin(¢)/,,.9 + @, cos(¢)cos(y)/,.¢ + w, sin(¢)sin(9)sin(y)/,.¢ 
~ 3a, sin(@)cos(9)sin(g)/_, + 3a? sin(@)cos(@) sin(¢)/,,, + h, 
+ ph, —y sin(@)h, —@, cos(@) sin(y )h, — Acos(¢)h, ~ycos(8)sin(¢)h, 
+@, cos(¢)cos(y )h, + @, sin(¢)sin(A) sin(y)h, — M,,) 


Then, after developing the simulink model, I analyzed the dynamics of the 
spacecraft as a function of time when a disturbance torque was applied. This can be seen 
in Figures 7-13 through 7-15; also, see Appendix G. Note that without control the 


spacecraft will become unstable very quickly. 
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Figure 7-13: Disturbance Torques Applied to the Spacecraft - Large Angle Scenario. 
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Figure 7-14: The three body components of angular velocity due to a disturbance torque 


with no control - Large Angle Scenario. 
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Figure 7-15: Spacecraft Angular Response Due to Disturbance Torque - Large 
Angle Scenario. 
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D. ADCS CONTROL AND PERFORMANCE ANALYSIS - YAW 
STEERING 


In order to keep the antenna beam pointing accuracy and the solar array panel 
pointing accuracy between acceptable limits (0.1° and 1° respectively), a control system 
must determine the required movements of the reaction wheel. This controller (control 
loop) could be developed in many different ways. In our case the angle and the angular 


speed measured by the sensors are used. In this way the system corrects the deviation 


based on the actual deviation and the rotation speed. The control loop is given in Figure 


7-16. 


0,0, y/ 


oO Satellite 
P rer oO per > W per 


PID 
Controller 









Figure 7-16: The Control Loop 


The block scheme of the system model is shown in Figure 7-17. The satellite 
dynamics are described by the equations 7-16. The BAPTA controller controls the actions 
of the BAPTA, which again influences the motions about the y-axis. There are three RW 
controllers each controlling a reaction wheel. The OBC has two sensor inputs, the orbital 


angle fi and the angle between the orbital plane and the Sun-Earth-vector 6. 
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Figure 7-17: System Model 


For the computer simulation Matlab and Simulink were used (see Appendix G). 
And, three situations were considered. First, a step input due to thruster misalignment. 
Second, the sinusoidal input due to solar pressure torques; and third, an impulse input of 
approximately one order of magnitude larger than the others (see simulink diagrams and 
Matlab code in Appendix G). Several combinations of the three scenarios were also 
tested. When modeling the simulations, disturbance torque values and scenarios were 
used that rarely would occur on the satellite. This was done in an effort to show the 
robustness of the satellite control. In general, we are trying to look at the worst case 


scenarios. Also, the control values selected are indicated in Table 7-2. 
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Table 7-2: Control Law Values 


The first simulation imparted impulse torques along the roll, pitch and yaw axes. 
The magnitudes selected were greater than those expected to actually occur on orbit and it 
be would extremely unlikely to have 3 simultaneous impulse torques. As can be seen in 
Figure 7-18, the impulse torques occurred at time, t=100 seconds. Figure 7-20 illustrates 
the angular responses for roll, pitch, and yaw. A thorough study of the responses reveal 
the greatest error occurs from the effects of yaw steering and not the impulse torques. 
And, the required pointing accuracies have been maintained throughout the orbit. Next, 
with a step input along all three axes, the angular responses were evaluated and found to 
be within operational limitations. It should be noted that the greatest error occurs during 
this scenario. This is because the step disturbance is occurring at the same time the yaw 
inversion is taking place. So in effect, we have an additive disturbance among the two. 
Then, a sinusoidal disturbance torque was applied along all three axis simultaneously. 
This can be seen in Figures 7-21 through 7-23. And, once again the pointing accuracies 
fell within limits. Finally, a simulation of sinusoidal, impulse, and step disturbance torques 
were applied along the roll, pitch, and yaw axes respectively. And, in this scenario all 
pointing accuracy requirements were maintained. 

Throughout the simulations the maximum error occurred when the yaw inversion 
took place. It is at this point that the yaw angular velocity is greatest. This is due to the 
fact that the yaw angle and angular rate are coupled into the roll and pitch. So, when yaw 
rate of change is greatest we see our greatest error along all three attitude axes. Also, the 
yaw angular velocity and angular acceleration is not continuous at their maximum and 
minimum values. This in turn creates a rapid change condition which then allows for more 
error to be introduced. 

The final simulation was a sinusoidal input along the pitch axis and was performed 


to illustrate the value of disturbance torque that would cause the system to exceed 
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pointing accuracy requirements. As can be seen in Figures 7-24 through 7-26, a sinusoidal 
input of 0.01 N-m caused the pitch axis pointing accuracy to exceed 0.1°. Recall that the 


maximum disturbance torque anticipated is 0.001 N-m. Therefore a built in safety margin 


exists. 
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Figure 7-18: Impulse Disturbance Torque for Control of Yaw Steering Scenario #1 


(Time in seconds). 
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Figure 7-19: Spacecraft Angular Velocities During the Control of Yaw Steering Scenario 


#1 (Time in seconds). 
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Figure 7-20: Spacecraft Angular Response Due to Impulse Disturbance Torque - Yaw 


Steering Scenario #1 (Time in seconds). 
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Figure 7-21: Sinusoidal Disturbance Torque for Yaw Steering Scenario #2. 


x 10° Angular Velocities vs Time 


phidot (deg/sec) 






thetadot (deg/sec) 


psidot (deg/sec) 


Time (sec) 4 


Figure 7-22: Spacecraft Angular Velocities During the Control of Yaw Steering Scenario 
#2. 
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Figure 7-23: Spacecraft Angular Response Due to Sinusoidal Disturbance Torque - Yaw 


Steering Scenario #2. 
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Figure 7-24: Sinusoidal Disturbance Torque for Yaw Steering Scenario #3. 
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Figure 7-25: Spacecraft Angular Velocities During the Control of Yaw Steering Scenario 
#3. 
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Figure 7-26: Spacecraft Angular Response Due to Sinusoidal Disturbance Torque - Yaw 


Steering Scenario #3. 


135 


This Page Intentionally Left Blank 


136 





Vill. CONCLUSION 


A. SUMMARY OF DESIGN PROCESS 


In this thesis the preliminary design of an Attitude Dynamics and Control 


Subsystem for a Medium Earth Orbit communications satellite was performed. The 


process focused on the six critical steps in designing an ADCS which are stated as follows: 


N A S&S W LW 


. Design control modes. 

. Define or derive system level requirements. 

. Select type of spacecraft control by attitude control mode. 
. Quantify the disturbance environment. 

. Select and size ADCS hardware. 


. Define determination and control algorithms. 


Also presented was the important aspect of how the ADCS is inter-related with the 


other spacecraft subsystems. And, how many of the requirements and design drivers of 


the ADCS originate in the other spacecraft subsystems. 


The highlights and conclusions of the ADCS design process are summarized as 


follows. First, the control modes of the spacecraft were determined and defined. The 


control modes determined to be most appropriate for our satellite were the following: 


Co.UmhMNL)])LUCON eee lOCWsaliKND“D 


. Orbit insertion mode 

. Acquisition mode 

. Sun acquisition/sun hold mode 
. Earth acquisition mode 

. Normal mode 

. Slew mode 

. Contingency mode 


. Orbit transition mode 
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The second design milestone was the determination of system level requirements. 


Ten ADCS requirements were presented which are summarized below as follows: 


1. The ADCS must be fully operational without the aid of ground base control. 

2. The operational lifetime of the ADCS must meet or exceed ten years. 

3. The ADCS must maintain a UHF antenna beam pointing accuracy of less than 
or equal to 0.1°. 

4. The ADCS must have the capability to provide continuous yaw steering and 
solar array tracking such that the solar array pointing accuracy is better than or equal to 
a 

5S. The spacecraft must be fully operational during eclipse. 

6. The spacecraft must be stabilized and under control during all phases of the 
mission. 

7. The ADCS system should have system redundancy such that there are no single 
point failures. 

8. The ADCS should have a ground override capability. 

9. The spacecraft must have a cross-link antenna pointing accuracy of 0.07° which 
requires a spacecraft attitude knowledge requirement of 0.02° in roll, pitch, and yaw. 


10. The ADCS must have a positional knowledge requirement of + 350 meters. 


The next design milestone was determining the type of spacecraft control by 
attitude control mode and presenting a functional description of the ADCS. During this 
stage we determined what kind of sensors and actuators were to be used and then showed 
how they were to be used in each mode of operation. Also introduced in this stage of the 
design was the motion of the satellite and solar arrays due to the yaw steering 
requirement. 

Then, the design focused on quantifying the disturbance environment. In this 
process the derivation of the Moments due to solar pressure, gravity gradient, 
aerodynamic and geomagnetic were presented. And, an estimated disturbance value was 


attained. Next, the selection of the ADCS hardware was performed. Each sensor and 
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actuator was described and the reasons for selection of each component were given. The 
sensor and actuator selection was driven solely by the requirements initially stated. 

The final step in the design process was to define determination and control 
algorithms. Chapter VI examined yaw controlling and presented explanations dealing with 
control strategy and control software. Then, in chapter VII several simulations of the 


ADCS were presented and the development of the control law was explained. 


B. SUMMARY OF SPACECRAFT ADCS DESIGN 


The ADCS chosen for the MUS satellite was a Zero Momentum three axis system 
which can provide any combination of slew maneuvers in the roll, pitch, and yaw 
directions. The ADCS can maintain Nadir pointing accuracies to within 0.1°. At the same 
time it can continuously rotate in the yaw direction while the solar arrays are 
simultaneously rotating in order to keep a perpendicular orientation towards the sun. This 
will greatly decrease the complexity of the solar array control mechanisms and reduce 
spacecraft mass. 

As previously mentioned the spacecraft ADCS has eight control modes. Each 
control mode has peculiar operating characteristics which enable the spacecraft to 
maintain operations for ten or more years. The operational requirements of the control 
modes were satisfied by the utilization of various ADCS components as presented im 


Figure 8-1. 
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Figure 8-1: Block Diagram of the Attitude Dynamics and Control Subsystem 


The ADCS components selected were those that satisfied the subsystem 


requirements as stated in the previous section. The components selected for the MUS 


spacecraft are listed in Table 8-1. 
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ADCS Control 
Processor Soke 25 W 


2 


. 
a 
Unit Z 3.45 kg 34 W 







Sun Sensor 


GPS Antenna 
Earth/Sun Sensor 


Table 8-1: Spacecraft ADCS Components 


After the ADCS component selection had taken place, the satellite bus design was 
able to proceed by incorporating the equipment. The final ADCS configuration for the 
MUS design is shown in Figure 8-2. 
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Figure 8-2: ADCS Component Layout 


Finally, a simulation of the spacecraft control was performed. When expected 
disturbance torques and yaw steering requirements were imposed on the model; pointing 
accuracy requirements were maintained. The limits of the subsystem were tested and it 
was concluded that it would require disturbance torques approximately one order of 
magnitude larger than what was expected to cause spacecraft operating limits to be 


exceeded. 


142 





ADCS 
AFSCN 
BAPTA 
CG 
CMG 
CPU 
COTS 
dB 
DoD 
DT 
DSCS 
DuSD - Space 
EHF 
EPS 

ES 
ES/SSA 
GEO 
GG 
GHz 
GPS 
HF 
IRU 
Kbps 
LEO 
MEO 
MHz 
MILSATCOM 
MOI 
MSS 
MUS 
Nm 
Nms 
PMOI 
Pn 

RF 

RW 

SA 
SCOC 
Sec 
SHF 
SOCC 


APPENDIX A 


ACRONYMS 


Attitude Dynamics and Control System 
Arr Force Satellite Control Network 
Bearing and Power Transfer Assembly 
Center of Gravity 

Control Moment Gyro 

Central Processing Unit 

Commercial Off The Shelf 

Decibels 

Department of Defense 

Disturbance Torque 

Defense Satellite Communications System 
Deputy Undersecretary of Defense for Space 
Extremely High Frequency 
Electrical Power Subsystem 

Earth Sensor 

Earth Sensor / Sun Sensor Assembly 
Geostationary Earth Orbit 

Gravity Gradient 

Gigahertz 

Global Positioning System 

High Frequency 

Inertial Reference Unit 

Kilobits per second 

Low Earth Orbit 

Medium Earth Orbit 

Megahertz 

Military Satellite Communication 
Moments of Inertia 

Mobile Satellite Services 

Mobile Users Subscribers 

Newton Meter 

Newton Meter Second 

Principal Moments of Inertia 
Primary 

Radio Frequency 

Reaction Wheel 

Solar Arrays 

Satellite Control Operations Center 
Secondary 

Super High Frequency 

Satellite Operational Control Center 
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SS Sun Sensor/Star Sensor 


TT&C Telemetry, Tracking and Control 
UFO UFH Follow On 

UHF Ultra High Frequency 

VHF Very High Frequency 
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APPENDIX B 
GLOSSARY 


Angular Momentum. The vector product of the position vector and linear velocity of a 
particle in motion relative to an axis. 


Cyclic Torque. Torque that repeats itself in cycles such that a build up of angular 
momentum does not occur. 


Desaturation, The process of removing built up angular momentum from the reaction 
wheels. 


Disturbance Torque. A torque applied to the spacecraft such that its attitude is 
disturbed. 


Geostationary Orbit. A geosynchronous orbit with an inclination of 0° or near 0°, 
positioning the satellite directly above the equator. [Ref. 16: p. 2] 


Geosynchronous Orbit. An orbit with a period equal to the earth's rotational period of 
1,436 minutes and an altitude above the earth of ~ 19300 nm. Most geosynchronous 
satellites move in near circular, near equatorial orbits. [Ref. 16: p. 2] 


Global Coverage. Coverage of all latitudes and longitudes, 90° N to 90° S. 


ICO. A medium earth orbit constellation consisting of 10 active satellites and 2 on orbit 
spares, which will operate at 10,400 km, providing worldwide coverage. The system will 
consist of 2 orbital planes with an inclination of 45°. Proposed by ICO Global 
Communications, IOC for the system is planned for the 4th quarter of 1999, with FOC 
scheduled for the 3rd quarter of 2000. ICO, which will use TDMA for modulation, will 
support 4,500 simultaneous handheld users per satellite, and can also support aeronautical, 
fixed, mobile and pager users. Satellite service life is projected at 10 years. The system 
will operate in S-band and C-band, supporting data rates of 2.4 or 4.8 kbps. Space 
segment cost is estimated to be 2.6 Billion (including 7 earth stations and 3 gateway 
switches), estimated handset cost is $1,000, and projected operating cost is $1.00- 
$2.00/minute in addition to a $40.00 per month access charge. [Ref. 16: p. 3] 


Low Data Rate. For the purposes of the Mobile User Study, 64 Kbps and below. 
[Ref. 16: p. 4] 
Medium Earth Orbit. Implies an orbital altitude between Low Earth Orbit (LEO ~ 


100nm to 540nm) and Geosynchronous Orbit (GEO ~ 18225nm). Orbit of choice for 
many proposed Mobile Satellite Systems. [Ref. 16: p. 5] 
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Mobile Satellite Service (MSS). The military use of satellite based commercial 


communications service providing digital voice, data, paging or fax services to users 
equipped with small (typically handheld) terminals. For the purpose of the MUS, this 
definition refers to Military use of PCS. 

[Ref. 16: p. 5] 


Mobile User. A highly mobile and or disadvantaged user that requires a small terminal 
and or antenna (includes hand held terminals and terminals for platforms such as trucks, 
submarines, small ships, boats, and aircraft). [Ref. 16: p. 5] 


Nadir. An earth pointing orientation or in the direction of the earth. 


Odyssey. A medium earth orbit constellation consisting of 12 active satellites with spares 
launched as required, which will operate at 10,355 km (circular orbit), providing 
worldwide coverage at all latitudes. The system will consist of 3 orbital planes with an 
inclination of 50°. Proposed by TRW, the system is planned for IOC in the 4th quarter of 
2001, with FOC planned for the 4th quarter of 2002. ODYSSEY will use CDMA for 
modulation, and will support 9,500 simultaneous handheld users per satellite, as well as 
fixed and transportable users. Space segment cost 1s projected to be 3.2 Billion, estimated 
handset cost is $500-$700, and projected operating cost is less than $1.00/minute. [Ref. 
16: p. 6] 


Personal Communications System (PCS). Satellite based commercial communications 
service providing digital voice, data, paging or fax services to users equipped with small 


(typically handheld) terminals. [Ref. 16: p. 6] 


Polar Coverage. Coverage that extends to the extreme northern and southern polar 
regions (90° N to 90° S). [Ref. 16: p. 6] 


Secular Torque. A torque that is continuous thus causing a build up of Angular 
Momentum to occur. 


Slew. A reorientation of the spacecraft about the roll, pitch or yaw axes. 


Bl - Area. 


f - The angle in the orbital plane between the point with minimum distance (sub solar 


point) to the sun and the satellite. 


B - The Geomagnetic Field. 


C, - Coefficient of drag. 
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d, - Vectors from the center of mass to each external piece's center of pressure. 


M , - Moment due to solar radiation pressure. 


ae 


p, - Photons that are specularly reflected from a surface. 





P, ~ Photons that are diffusely reflected from a surface. 
P, ~ Photons absorbed by a surface. 
p - Atmospheric density. 


oO - The angle between the orbital plane and the direction of the sun. 


0, - Critical value of 6 where an alternative control scheme must take place. 
F- - Radiation force due to absorbed photons. 
F’, - Force due to solar radiation pressure. 


BF - Radiation force due to specularly reflected photons. 





i - Satellite orbit inclination. 


I, - Vector distance from center of mass to center of pressure. 


n - Solar Array angle of rotation. 


u,, - Unit vector in velocity direction. 





s - Solar radiation unit vector. 


Q.- Right Ascension of the ascending node. 
T - Spacecraft torque due to areodynamic effects. 
T® - Spacecraft torque due to gravity effects. 


T’” - Spacecraft torque due to the magnetic field effects. 


147 


T.. - Torque about the spacecraft's roll axis. 


es 


4 - Torque about the spacecraft's pitch axis. 


T. - Torque about the spacecraft's yaw axis. 
v - Earth centered true anomaly. 


X. - The satellite's roll axis. 


Y. - The satellite's pitch axis. 


Z, — axis - Yaw axis or axis parallel to spacecraft UHF antenna 


y - Yaw rotation angle. 
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APPENDIX C 


MATLAB CODE FOR YAW STEERING 


% Danny Busch 21 August 1997 
% AA 4871 

% This .m file plots the Yaw Rotation Angle and Solar 

% Array Rotation angle as a function of time and as 

% a function of the Sun Rotation Angle in order to 

% maintain solar array tracking. 


a=17378; 
u=398600.5; 
t=0:1:22798; 
x=(u/(a%3))*(1/2); 
B=(x)*t; 


figure(1) 

d=1*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/pi); 
subplot(2, 1,1); 
plot((B*(180/pi)),Psi); 


hold on 

d=5*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/p1); 
plot((B*(180/pi)),Psi); 


d=15*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/p1); 
plot((B*(180/pi)),Psi); 


d=30*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/p1); 
plot((B*(180/pi)),Ps1); 


d=45*(p1/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/pi); 
plot((B*(180/pi)),Psi); 
d=60*(p1/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/pi); 
plot((B*(180/pi)), Psi); 


d=75*(pi/180); 
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Psi=(atan((-sin(B )/tan(d))))*(180/pi1); 
plot((B*(180/pi)),Psi); 


d=90*(pi/180); 

Psi=(atan((-sin(B )/tan(d))))*(180/pi); 
plot((B*(180/pi)),Psi), grid; 

xlabel(‘Sun Rotation Angle (degrees)’); 
ylabel(‘Yaw Command Angle (degrees)’); 
title("Nominal Yaw Trajectories’); 
axis([0,360,-100,100]); 


gtext(‘d=1’); 
gtext('d=5'); 
gtext(‘'d=15'); 
gtext('d=30"); 
gtext('d=45'); 
gtext('d=60'); 
gtext('d=75'); 
gtext('d=90'); 
hold off 


subplot(2, 1,2); 

d=1*(pi/180), 

Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/p1)),(n*(180/pi))); 

ylabel(‘Solar Array Rotation Angle in degrees’); 
xlabel('Sun Rotation Angle in Degrees’); 
title(Solar Array Trajectories’); 

hold on 


d=5*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))); 


d=15*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(1 80/pi))); 


=30*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))); 
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d=45*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(1 80/p1)),(n*(180/p1))); 


d=60*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Ps1)./tan(B))); 
plot((B*(180/pi)),(n*(180/p1))); 


d=75*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))); 


d=90*(pi/180): 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot((B*(180/pi)),(n*(180/pi))), grid; 
axis({0,360,-100,100]); 


gtext(‘d=1'); 

gtext('d=5'); 

gtext('d=15'); 
gtext('d=30'); 
gtext('d=45'); 
gtext('d=60'); 
gtext('d=75'); 
gtext('d=90'); 


figure(2) 
t=t*(1/60); % Convert time to minutes 


d=1*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/p1); 
subplot(2,1,1); 

plot(t,Psi); 


hold on 

d=5*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/p1); 
plot(t,Psi); 


d=15*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/pi); 
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plot(t,Psi); 


d=30*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/pi); 
plot(t,Psi); 


d=45*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/p1); 
plot(t,Psi); 


d=60*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/pi); 
plot(t,Psi); 


d=75*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/p1); 
plot(t,Psi); 

axis([{0,380,-100,100]); 


d=90*(pi/180); 
Psi=(atan((-sin(B)/tan(d))))*(1 80/pi); 
plot((B*(180/p1)),Psi), grid 

xlabel("Time in minutes’); 

ylabel("Yaw Command Angle (degrees)'); 
title(‘Nominal Yaw Trajectories’); 


Yogtext(’d=1'); 
Ygtext('d=5'); 
Ygtext('d=15'); 
Yogtext('d=30'); 
Ygtext('d=45'); 
Yogtext('d=60'); 
Yogtext('d=75'); 
“ogtext('d=90'); 
hold off 


subplot(2, 1,2); 

d=1*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))); 

ylabel(‘Solar Array Rotation Angle in degrees’); 
xlabel(‘Time in minutes’); 

title(Solar Array Trajectories’); 

hold on 
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d=5*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))), grid; 


d=15*(pi/180), 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/p1))); 


d=30*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))); 


d=45*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))); 


d=60*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Ps1)./tan(B))); 
plot(t,(n*(180/pi))); 


d=75*(pi/180); 
Psi=(atan((-sin(B )/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/pi))); 


d=90*(pi/180); 
Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Psi)./tan(B))); 
plot(t,(n*(180/p1))); 
axis({0,380,-100,100}); 


Ygtext('d=1'); 

Ygtext('d=5'); 

“gtext('d=15'); 
Ygtext('d=30'); 
“Ygtext('d=45'); 
Yogtext('d=60'); 
Yogtext('d=75'); 
“gtext('d=90'); 
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figure(3) 


subplot(2,1,1); 

d=5*(p1/180); 
Psi=(atan((-sin(B)/tan(d))))*(180/pi); 
plot(t,Psi); 


for 1=1:length(Ps1); 
Psi(i)=90; 
end; 


hold on 

plot(t,Psi), grid; 

xlabel('Time in minutes’); 

ylabel('Yaw Command Angle (degrees)'); 
title(‘Nominal Yaw Trajectories’); 

“%ptext('d=5'); 

Ygtext("Yaw Angle = Constant 90 degrees for d < 5’); 
axis({0,380,-140,140]) 


y=[120,-120,120,-120]; 
g=[t(1),t(10695),t(12045),t(22799)]; 


subplot(2, 1,2); 
plot(q,y) 


hold on 

d=5*(pi/180); 

Psi=(atan((-sin(B)/tan(d)))); 
n=atan(-(sin(Ps1)./tan(B))); 

plot(t,(n*(180/pi))), grid; 
axis({0,380,-150,150]); 

xlabel(‘Time in minutes’); 

ylabel(‘Solar Array Rotation Angle in degrees’); 
title(‘Solar Array Trajectory for d < 5 degrees’); 
Yogtext(‘'d=5'); 
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APPENDIX D 


MATLAB CODE TO COMPUTE DIRECTION COSINE MATRIX 


% Danny K. Busch 4 October 1997 
% Thesis 

% 

% This is a supplemental .m file to compute the initial dem 

% from the orbit/nominal frame to the principal axis frame. 

% I am making the assumption that the orbit/nominal frame 

% is lined up body axes when this starts. So this is the initial 
% dcm from both the body axes and orbit frame to the principal 
% axes. 


% From the spreadsheet we have the Inertia Matrix. 
Jold=[13337 -28.97 0;-28.97 14083 17.72;0 17.72 14113]; 
[evec, princ]=eig(Told) 


% Where princ is the matrix containing the principal moments 
% of inertias, and Evec is the transpose of the dcm to get there. 


dcm=evec' 
% To prove this, multiply dcm*Iold*dcm' to get from the old 
% (non-principal) system to the new (principal) system. 


Inew=dcem*Iold*dcm' 
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% Danny K. Busch 4 October 1997 
% Thesis 


% Plot the body frame vs the Principal Inertia Axis 


dcm=[0.9992,0.0388,0.0009;0.0354,-0.9026,0.4290;-0.0158,0.4287,0.9033] 
x=[0,0,0;dcm(1,1),dem(1,2),dcem(1 ,3)]; 
y=[0,0,0;dmc(2,1),dem(2,2),dcm(2,3)]; 
z=[0,0,0; dcm(3,1),dcm(3,2),dem(3,3)]; 


u=[0,0,0; 1,0,0]; 
v=[0,0,0;0, 1,0]; 
w=[0,0,0;0,0, 1]; 


figure(1); 

plot3(x(:,1),x¢.,2),x(:,3),'r5yC, 1), vC,2), vC,3),'r,Z€, 1),z¢,2),ZC,3),'r) 
hold on 

plot3(u(:,1),u(.,2),uC,3),’b5 VC, 1),vC,2),vC,3),'bw, 1),wl,2),w,3),'b') 
hold off 

grid; 

xlabel('X’), ylabel("Y'),zlabel('Z’); 

title(‘Body Frame vs Principal Inertial Frame’); 
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APPENDIX E 


SIMULATION OF THE TOTAL EXTERNAL DISTURBANCES 


% Danny K. Busch 27 July 1997 

% Thesis 

% ADCS Design 

% 

% Gravity Gradient + Solar Pressure + Aerodynamic Pressure Torques 
% 

% Simulates the attitude motion of a s/c vs time. 

% That is the orientation of the b frame wrt to the 

% orbit/nominal frame. Without attitude control. 

% 


function xdot=total(x); 
°% Define constants 
11=45.01961; 
12=35.53275: 
13=59.65711; 
Ky=(13-12)/I1; 
Kr=(13-11)/12; 
Kp=(12-11 )/13; 


% Semi-major axis of sun/earth system 
asun=1.49598e8; 


°% Sun constant 
us=1.327el11; 
to=0; 

Mos=pi; 


% Arguement of perigee 
w=0.4363323; 


% Longitude of the ascending node 
omega=(pi/2); 


% Inclination 
inc=(p1/4); 


% Eccentricity of Sun/Earth system 
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es=.016708; 


% Eccentricity of Earth/Satellite system 
ee=0; 


ue=3.986012e5; 


% Radius of the Earth 
Re=6378; 


% Generaic cause normalize later, vector from sun to earth 
r=10; 


°% R=dist fm center of earth to center of mass of s/c 
% User friendly conversions 


w1l=x(1); 
w2=x(2); 
w3=x(3); 
all=x(4); 
al2=x(5); 
al3=x(6); 

a2 1=x(7); 
a22=x(8); 
a23=x(9); 

a3 1=x(10); 
a32=x(11); 
a33=x(12); 
rl=x(13); 
r2=x(14); 
r3=x(15); 
rdot 1=x(16); 
rdot2=x(17); 
rdot3=x(18); 
time=x(19); 


R=sqrt(r142+12/2+132); 
v=sqrt(rdot1“2+rdot2“2+rdot32); 
nudot=v/R; 


Yo[eRs]|h=r vector fm sun to earth expressed in heliocentric 
Yotrame 
“ol[eRs|I=r vector fm sun to earth expressed in earth 
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“% centered inertial 

“e[eRs]o= r vector fm sun to earth expressed in the 
“orbit frame 

%[eRs]b= r vector fm sun to earth expressed in the 
Ybody frame 


% Find nu and d as a function of time 


Moearth=0; 

toearth=0; 

Yoa = semimayjor axis of satellite earth system 

a=2 1000; 

ne=(ue/a%3)(1/2); 

% ne mean orbital rate 

%M-Mo=n(t-to) general equation mean anomoly 


Me=ne"* time; 

Yonu earth 

nu=Me + 2*ee*sin(Me) + 1.25*((ee)“2)*sin(2*Me); 
ns=sqrt(us/asun’3); 

“oMs-Mos=ns(time-to); mean anomoly for earth sun system; 
Ms=ns*time+p1; 

Yoorbital parameter nu for sun earth system 
d=Ms+2*es*sin(us)+1.25*(es*2)*sin(2*Ms); 


%Now we have nu and d_ or nu for sun and earth 
% R vector from sun to earth in earth centered inertial 


eRsI=[-r*cos(d);-r*cos(.4084)*sin(d);-r* sin(.4084)* sin(d)]; 
% multiply by -1 so we R points from earth to sun according to Agrwal 


% Define dem from earth centered inertial to rotating orbital frame 
J = [cos(nu) sin(nu) 0;-sin(nu) cos(nu) 0;0 0 1); 

K = [cos(w) sin(w) 0;-sin(w) cos(w) 0;0 0 1]; 

L =[1 0 0;0 cos(inc) sin(inc);0 -sin(inc) cos(inc)]; 

Q = [cos(omega) sin(omega) 0;-sin(omega) cos(omega) 0;0 0 1]; 
oCi = J*FK*L*Q; 


% Now transform to the orbit frame 


eRso = oCi*eRsI; 
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%define dcm from orbit to body frame 

bCo = [al1 a21 a31;a12 a22 a32;a13 a23 a33]; 
% Now tranform R vector to orbit frame 

eRsb = bCo*eRso; 


% Now compute the unit vector from earth to sun 
Shatb = eRsb/(norm(eRsb)); 
s= Shatb; 


% Now we have S hat expressed in the b frame as a 
% function of d and nu, which is a function of time 


% Define n hat vectors in the body frame 


nl = [0; 1; 0}; 

n2 = [0.866; 0.5; 0]; 
n3 = [0.866; -0.5; 0]; 
n4 = [0; -1; 0]; 

n5 = [-0.866; -0.5; 0]; 
n6 = [-0.866; 0.5; 0]; 
n7 = [0; 1; 0]; 

n8 = [0; -1; 0]; 

n9 = [-1; 0; 0]; 

nl0 = [0; 0; 1}; 

nl] = [0; 0; -1]; 

nl2 = [0; 0; -1]; 

n13 = [0; 0; 1}; 


% Now define the dr vectors from center of mass to center 
% of pressure 


dr1 = [0.11; 0.5; 0]; 
dr2 = [0.543; .25; 0}; 
dr3 = [0.543; -.25; 0]; 
dr4 = [0.11; -0.5; 0]; 
dr5 = [-0.543; -0.25; 0]; 
dr6 = [-0.543; 0.25; 0]; 
dr7 = [-0.756; 0.5; 0]; 
dr8 = [-0.756; -0.5; 0]; 
dr9 = [-0.39; 0; 0]; 
dr10 = [-.756; 0; 0.5]; 
dr11 = [-.756; 0; -0.5]; 
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aileze Ot 1: 0; -0.5]; 
dr13 = [0.11; 0; 0.5]; 


Fe = 1358; 

C = 2.9979E8; 
P = Fe/C; 

A = zeros(1,13); 


% F = PA(nhat dot shat)(2*nhat) 
% Calculate the (nhat dot shat) 


hat = zeros(1,13); 
hat(1) = sum(n1.*s); 
if hat(1)<0; 
hat(1)=0; 
else 
hat(1)=hat(1); 
end 


hat(2) = sum(n2.*s); 
if hat(2)<0; 
hat(2)=0; 
else 
hat(2)=hat(2); 
end 


hat(3) = sum(n3.*s); 
if hat(3)<0; 
hat(3)=0; 
else 
hat(3)=hat(3); 
end 


hat(4) = sum(n4.*s); 
if hat(4)<0; 
hat(4)=0; 
else 
hat(4)=hat(4); 
end 


hat(5) = sum(n5.*s); 
if hat(5)<0; 
hat(5)=0; 
else 
hat(5)=hat(5); 
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end 


hat(6) = sum(n6.*s); 
if hat(6)<0; 
hat(6)=0; 
else 
hat(6)=hat(6); 
end 


hat(7) = sum(n7.*s); 
if hat(7)<0; 
hat(7)=0; 
else 
hat(7)=hat(7); 
end 


hat(8) = sum(n8.*s); 
if hat(8)<0; 
hat(8)=0; 
else 
hat(8)=hat(8); 
end 


hat(9) = sum(n9.*s); 
if hat(9)<0; 
hat(9)=0; 
else 
hat(9)=hat(9); 
end 


hat(10) = sum(n10.*s); 
if hat(10)<0; 
hat(10)=0; 
else 
hat(10)=hat(10); 
end 


hat(11) = sum(n11.*s); 
if hat(11)<0; 
hat(11)=0; 
else 
hat(] 1)=hat(1 1); 
end 


hat(12) = sum(n12.*s); 
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if hat(12)<0; 
hat(12)=0; 

else 
hat(12)=hat(12); 

end 

hat(13) = sum(n13.*s); 

if hat(13)<0; 
hat(13)=0; 

else 
hat(13)=hat(13); 

end 


% hat = nhat(i) dot shat(1) 
% hat = [-0---00--0--o] 


% Define the area matrix 


een >>). 331 331 1.77 331 33) .724 724% 
ps = .4; 
pd =.4; 


Yo From Agrwal’s Book p. 134 
% Ms = PA(n dot s)*r X [(1-ps)*s + 2(ps + 1/3pd)n] 


Ms = zeros(3, 13); 
“% Assign intermediate variables 


G = zeros(3, 13); 

G(:,1) = P*A(1)*hat(1)*dr1; 
G(:,2) = P*A(2)*hat(2)*dr2; 
G(.,3) = P*A(3)*hat(3)*dr3; 
G(:,4) = P*A(4)*hat(4)*dr4; 
G(:,5) = P*A(5)*hat(5)*dr5; 
G(:,6) = P*A(6)*hat(6)*dr6; 
G(:,7) = P*A(7)*hat(7)*dr7; 
G(.,8) = P* A(8)*hat(8)*dr8; 
G(:,9) = P*A(9)*hat(9)*dr9; 
G(:,10) = P*A(10)*hat(10)*dr10; 
G(.,11) = P*A(11)*hat(11)*dr11; 
G(:,12) = P*A(12)*hat(12)*dr12; 
G(.,13) = P*A(13)*hat(13)*dr13; 


T = zeros(3,13); 
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T(:,1) = ((1-ps)*s) + ((2*(ps + 0.333 *pd))*n1); 
T(:,2) = ((1-ps)*s) + ((2*(ps + 0.333*pd))*n2); 
T(.,3) = ((1-ps)*s) + ((2*(ps + 0.333*pd))*n3); 
T(.,4) = ((1-ps)*s) + ((2*(ps + 0.333 *pd))*n4); 
T(:,5) = (C-ps)*s) + ((2*(ps + 0.333 *pd))*n5); 
T(:,6) = ((1-ps)*s) + ((2*(ps + 0.333 *pd))*n6); 
T(:,7) = ((1-ps)*s) + ((2*(ps + 0.333*pd))*n7); 
T(:,8) = ((1-ps)*s) + ((2*(ps + 0.333*pd))*n8); 
T(:.,9) = ((1-ps)*s) + ((2*(ps + 0.333 *pd))*n9); 
T(:,10) = ((1-ps)*s) + ((2*(ps + 0.333*pd))*n10); 
T(:,11) = ((1-ps)*s) + ((2*(ps + 0.333*pd))*n1 1); 
T(.,12) = ((1-ps)*s) + ((2*(ps + 0.333*pd))*n12); 
T(:,13) = ((1-ps)*s) + ((2*(ps + 0.333*pd))*n13); 


Z = cross(G,T); 


% We now have a 3 by 13 matrix containing torques in b1 direction for 
% thirteen areas in the first row, etc... 


% Solar pressure moments in b direction 
Msp1 = sum(Z(1,:)); 
Msp2 = sum(Z(2,:)); 
Msp3 = sum(Z(3,:)); 


% (1/2)*p*Cd* A*v*v*(uvhat X Icp) 
% Define p as a function of R 


if R < 6578; 
p=4e-10; 
elseif R < 6678; 
p=5e-11; 
elseif R < 6778; 
p=1.5e-11; 
elseif R < 6878; 
Pa c- ize 
elseif R < 6978 
p=2e-12 
else 
p=8e-13; 
end 


% Define Cd as a constant given on P.90 of handout. 


Cd=3.2; 


164 


’% Define Area constants in the b1, b2 and b3 directions. This is the total 
% cross sectional area in the b1 hat direction,...b2 hat direction and b3 
% hat direction. 


Ab1=1.77; 
Ab2=.884: 


? 


Ab3=1.61; 

Abody=[Ab1 Ab2 Ab3]'; 

% Now build the velocity vector for the I frame 
RDOTi=[rdot1 rdot2 rdot3]'; 

% Now convert velocity from I frame to B frame 
RDOTb=bCo*0Ci*RDOTI; 

“% Now find the unit vector of RDOTb, uvhat. 
uvhatb=RDOTb./norm(RDOTb); 

% Now compute area perpendicular to uvhat. 
AT=abs(Abody. *uvhatb); 

% Compute total area 

ATOTAL=sum(AT); 

% Now compute (uvhatb X Icp) 

Icp=zeros(3,13); 

Icp=[drl dr2 dr3 dr4 dr5 dr6 dr7 dr8 dr9 dr10 dr11 dr12 dr13}; 
Maero=(1/2)*p*Cd* ATOTAL*v*v*(cross(uvhatb,Icp)); 
Maerototal=[sum(Maero(1,:)) sum(Maero(2,:)) sum(Maero(3,:))]'; 
Msp 1=sum(Z(1,:)); 

Msp2=sum(Z(2,:)); 

Msp3=sum(Z(3,:)); 

Ml1gg=((3 *ue)/R%3)*(I3-12)*al2*al3; 


M2ge=((3*ue)/R‘3 )*(I1-13)*al3*al 1; 
M3 gg=((3 *ue)/R%3)*(12-11 )*al 1 *al 2; 
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xdot(1)= ((Msp 1+Maerototal(])+M1gg)/I1) - Ky*w2*w3; 
xdot(2)= ((Msp2+Maerototal(2)+M2gg)/12) - Kr*w3* wl; 
xdot(3)= ((Msp3+Maerototal(3)+M3gg)/I3) - Kp*wl1*w2; 
xdot(4)=a12*(w3-a33*nudot)-al3*(w2-a32*nudot); 
xdot(5)=al3*(w1-a31*nudot)-al 1*(w3-a33*nudot); 
xdot(6)=al 1*(w2-a32* nudot)-al2*(w1-a3 1 *nudot), 
xdot(7)=a22*(w3-a33* nudot)-a23*(w2-a32*nudot); 
xdot(8)=a23*(w1-a3 1*nudot)-a2 1*(w3-a33*nudot); 
xdot(9)=a21*(w2-a32*nudot)-a22*(w1l-a3 1 *nudot); 
xdot(10)=a32*(w3-a33*nudot)-a33*(w2-a32*nudot); 
xdot(1 1)=a33*(w1-a3 1 *nudot)-a3 1 *(w3-a33* nudot); 
xdot(12)=a3 1*(w2-a32*nudot)-a32*(w1-a3 1 *nudot); 
xdot(13)=x(16); 

xdot(14)=x(17); 

xdot(15)=x(18); 

xdot(16)=((-ue)/R%3)*r1; 

xdot(17)=((-ue)/R‘%3)* 12; 

xdot(18)=((-ue)/R%3)*1r3; 
xdot(19)=Msp1+Maerototal(1)+M1 gg; 
xdot(20)=Msp2+Maerototal(2)+M2gg; 
xdot(21)=Msp3+Maerototal(3)+M3 gg; 
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APPENDIX F 


SMALL ANGLE SIMULATION 


a 
Pitch Response Simulation (f= Lire 7 


Clock To Workspace2 
TQuitch |To Workspace3 



















Steplinput 
Dre Ss 
lyy s<+3°wo2* (Ixx4zz) | 
, rad to degrees 10 Workspace 
Sine Wave 


Transfer Fen 


PID Controller 


Integrator To Workspace? 
impulsea 


impulseb 


Impulse 


169 


% Danny K. Busch 4 October 1997 
% Thesis 
% This is the initial conditions used with the simulink. 


wo=2.75578e-4; 
Ixx= 13337; 
lyy=14083; 
Izz=14113; 
Kp=.3; 

Tp=433; 

Ky=. 1; 

Ty=751; 

Kr=.5; 
Tr=2060; 
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% Danny Busch 4 October 1997 
% Thesis 

% This file plots torque, angle (theta), and reaction 

°% wheel momentum vs time. 


figure(1) 


subplot(3,1,1); 

plot(time, TQpitch), grid; 
title(‘Pitch Torque vs Time’); 
xlabel('Time (sec)'); 
ylabel('Torque (N-m)'); 
orient tall 


subplot(3,1,2); 

plot(time, theta), grid; 

title(Pitch Angle vs Time’); 
xlabel('Time (sec)'); 

ylabel(‘Pitch Angular Response (deg)’); 
orient tall 


subplot(3,1,3) 

plot(time,hp), grid; 

title(‘Pitch Wheel Angular Momentum vs Time’); 
xlabel('Time (sec)'); 

ylabel(‘Angular Momentum (N-m-s)’); 


orient tall! 


17] 


S 
Roll Yaw Response Simulation 


TQroll -|To Workspace3 
Clock To Workspace2 









Steplinput 
a ; : 
ae x%.$2+4°WOr2 (lyytz CMe 
Sine Wave rad todegrees 10 Workspace 
impulsea 
impulseb 
Integrator To Workspace? 
Step Fen 
ie! 
MY 
, rad to degrees To Workspace4 
7. evel Transfer Font 
dc cornp 
Env Torque 






impulseat 


PID Controlier1 
impulsebt a 
Impulse 1 — 
Gamn2 backup system 


es Ene 


Gaint Integrator? To Workspace5 
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% Danny K. Busch 4 October 1997 
% Thesis 

% Plot .m file used for thesis to aid in plotting 

% the results from Roll and Yaw 

% Plots torque, angle and momentum vs time. 


figure(1) 

subplot(3,1,1); 

plot(time,TQroll), grid; 

title(‘Roll Torque vs Time’); 
xlabel('Time (sec)'); 

ylabel("Torque (N-m)’); 

orient tall 

subplot(3, 1,2); 

plot(time,phi), grid; 

title(‘Roll angle phi’); 

xlabel('Time (sec)’); 

ylabel(‘Roll Angular Response (deg)'); 
orient tall 

subplot(3, 1,3) 

plot(time,hr), grid; 

title’ Angular Momentum vs Time’); 
xlabel('Time (sec)'); 

ylabel(’Angular Momentum (N-m-s)’); 
orient tall 

figure(2) 

subplot(3,1,1); 

plot(time, TQyaw), grid; 

title("Yaw Torque vs Time’); 
xlabel('Time (sec)'); 

ylabel(‘Torque (N-m)'); 

orient tall 

subplot(3,1,2); 

plot(time,psi), grid; 

title("Yaw angle ps1’); 

xlabel('Time (sec)'); 

ylabel("‘Yaw Angular Response (deg)’); 
orient tall 

subplot(3, 1,3) 

plot(time,hy), grid; 

title’ Angular Momentum vs Time’); 
xlabel('Time (sec)'); 

ylabel('‘Angular Momentum (N-m-s)’); 
orient tall 
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APPENDIX G 


LARGE ANGLE SIMULATION 


“% Danny K. Busch 3 October 1997 
% Thesis 

% 

7% Using S Function to model the dynamics of the Spacecraft 


% Define S-Function 
function[sys, xo]=large(t,x,u,flag) 
% VDPM EOM as an M-file 


Ixx=13337; 
lyy=14083; 
Izz=14113; 
wo=2.75578e-4; 


if abs(flag)==1, 
% Return rates 


sys(1)=x(2); 

sys(2)=(- 1/Ixx)*(x(4)*cos(x(1))*Izz*x(6)+x(6)*x(6)*cos(x(3))*sin(x(1))*Izz- 
wo*cos(x(1))*cos(x(5))*Izz*x(6)- 
wo*sin(x(1))*sin(x(3))*sin(x(5))*1zz*x(6)+x(4)*x(4)*sin(x(1))*Iyy- 
x(6)*cos(x(3))*cos(x(1))*lyy*x(4)- 
wo*sin(x(1))*cos(x(5))*lyy*x(4)+wo*cos(x(1))*sin(x(3))*sin(x(5))*lyy*x(4)+u(4)+x(4) 
*cos(x(1))*u(3 )+x(6)*cos(x(3))*sin(x(1))*u(3)-wo*cos(x(1))*cos(x(5))*uG )- 
wo*sin(x(1))*sin(x(3 ))* sin(x(5))*u(3 )+x(4)*sin(x(1))*u(2)- 
x(6)*cos(x(3))*cos(x(1 ))*u(2)- 
wo*sin(x(1))*cos(x(5))*u(2)+wo*cos(x(1))*sin(x(3))*sin(x(5))*u(2)- 
3*wo*wo*sin(x(1))*cos(x(1))*cos(x(3))*cos(x(3 ))*1zz+3 *wo*wo* sin(x(1))*cos(x(1))*c 
os(x(3))*cos(x(3))*Iyy); 


sys(3)=x(4); 

sys(4)=(-1/lyy)*(- 
x(4)*sin(x(1))*Ixx*x(2)+x(6)*cos(x(3))*cos(x(1 ))*Ixx*x(2)+wo*sin(x(] ))*cos(x(5))*Ixx 
*x(2)-wo*cos(x(1))*sin(x(3))*sin(x(5))*Ixx*x(2)- 
x(2)*1zz*x(6)+x(6)*x(6)*sin(x(3))*Izz+wo*cos(x(3))*sin(x(5))*Izz*x(6)- 
x(4)*sin(x(1))*u(1)+x(6)*cos(x(3))*cos(x(1))*uC1)+wo* sin(x(1))*cos(x(5))*u(1)- 
wo*cos(x(1))*sin(x(3 ))*sin(x(5))*uC1 )- 
x(2)*u(3)+x(6)* sin(x(3 ))*u(3)+wo*cos(x(3))*sin(x(5))*u(3)+u(5)- 
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3*wo*wo'*sin(x(3))* cos(x(3 ))*cos(x(1))*1zz+3 *wo* wot sin(x(3))*cos(x(3))*cos(x(1))*1 
XX); 


sys(5)=x(6); 

sys(6)=(-1/Izz)*(x(2)* lyy*x(4)-x(6)*sin(x(3))*Iyy*x(4)- 
wo*cos(x(3))*sin(x(5))*lyy*x(4)-x(4)*cos(x(1))*Ixx*x(2)- 
x(6)*cos(x(3))*sin(x(1))* Ixx*x(2)+wo*cos(x(1))*cos(x(5))*Ixx*x(2)+wo*sin(x(1))*sin( 
x(3))*sin(x(5))*Ixx*x(2)+x(2)*u(2)-x(6)* sin(x(3))*u(2)-wo*cos(x(3))* sin(x(5))*u(2)- 
x(4)*cos(x(1))*u(1)- 
x(6)*cos(x(3))*sin(x(1))*u(] )+wo*cos(x(1))*cos(x(5))*u(1)+wo* sin(x(1))*sin(x(3))* sin( 
x(5))*u(1)+u(6)- 
3*wo*wo*sin(x(3))*cos(x(3))*sin(x(1 ))*Ixx+3* wo*wo* sin(x(3))* cos(x(3))* sin(x(1))* ly 
y); 


elseif flag==0 
% Return initial conditions 
sys=[6;0;6;6;0;0]; 


% [continuous states, discrete states, outputs, inputs, discontinuous roots, 
% direct feed through] 


elseif flag==3 
% If flag = 3, return system outputs 


sys(1)=x(1); 
sys(2)=x(2), 
sys(3)=x(3); 
sys(4)=x(4); 
sys(5)=x(5); 
sys(6)=x(6); 


else 


sys=[]; 
end; 
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% Danny K. Busch 4 October 1997 
% Thesis 
% This is the initial conditions used with the simulink. 


wo=2.75578e-4; 
Ixx=13337; 
lyy=14083; 
Izz=14113; 


Kr1; 


Kp=2.0; 
Ky=.2; 


%Kr=0.5; 
Tr=2060; 


“oK p=.3; 
Tp=433; 


%Ky=.1; 
Ty=751,; 
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% Danny Busch 13 October 1997 
% Thesis 

% 

°% .m file used to simulate the motion of the satellite 

function act=psicmd(t) 
act=(-1.000507)*atan(11.43005*sin((0.0002755937)*(t))); 


end; 
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